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ABSTRACT  OF  THESIS 


DYNAMIC  ADJUSTMENT  OF  MESOSCALE  CONVECTIVE 
LOWER-STRATOSPHERIC  OUTFLOWS 


Recent  observational  studies  of  upper-tropospheric  and  lower-stratospheric  winds  atop 
mesoscale  convective  systems  show  the  development  of  anticyclonic  outflow.  We  propose 
that  the  anticyclone  formation  can  be  partially  explained  by  the  gradient  adjustment  that 
follows  the  nearly  instantaneous  vertical  redistribution  of  mass  by  the  convection.  With 
the  convection  idealized  by  an  impulsive,  diabatic  mass  transfer  from  the  lower  troposphere 
to  a  layer  near  the  tropopause,  the  adjustment  process  is  examined  using  an  invertibility 
principle  developed  from  the  quasi-static  primitive  equations  for  axisymmetric,  inviscid, 
adiabatic  flow  on  an  /-plane  in  potential  radius  and  entropy  coordinates.  The  invertibility 
principle  is  solved  as  a  single,  nonlinear,  elliptic  problem.  Solutions  show  the  development 
of  an  anticyclone  aloft  with  cold  and  warm  temperature  anomalies  above  and  below, 
respectively.  Sensitivity  studies  indicate  that  the  anticyclone  strength  is  greatest  for  lower- 
stratospheric  injections  at  high  latitudes  including  the  effects  of  cloud-top  cooling.  As 
the  magnitude  of  the  anticyclone  increases,  the  inertial  stability  of  the  system  is  reduced, 
resulting  in  a  decrea^sed  partitioning  of  the  initial  available  potential  energy  to  the  balanced 
state  of  the  system. 
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Chapter  1 


INTRODUCTION 

Upon  witnessing  the  immensity  of  a  mesoscale  convective  system  (MCS)  and  observ¬ 
ing  the  vastness  of  its  cirrus  anvil,  we  immediately  wonder  what  influence  such  intense 
disturbances  have  on  the  flow  aloft.  Ultimately,  we  must  question  what  errors  the  dis¬ 
turbances  introduce  into  our  model  forecasts,  since  most  large-scale  numerical  weather 
prediction  (NWP)  models  do  not  resolve  mesoscale  convective  fluxes  of  mass  and  momen¬ 
tum.  Emphasizing  this  point,  Maddox  (1980a)  explains  that  mesoscale  convective  com¬ 
plexes  (MCC)  in  particular  are  “not  explicitly  predicted  by  operational  numerical  models 
even  though  they  are  shown  to  be  organized  in  a  distinctly  nonrandom  mode  on  scales 
that  cannot  be  considered  subgrid.”  This  statement  certainly  applies  not  only  to  MCCs, 
but  to  squall  lines  and  other  forms  of  MCS  as  well.  Lilly  (1975)  illustrates  the  point,  stat¬ 
ing  that  mesoscale  convection  has  “a  degree  of  organization  which  is  clearly  nonrandom 
but  also  largely  unresolvable  from  conventional  data  processed  in  conventional  operational 
ways.  After  its  organization  the  convective  storm  array  develops  very  high-energy  density 
and  one  may  suspect  that  .  .  .  the  large-scale  environment  loses  control  and  the  tail 
wags  the  dog  (emphasis  added).”  No  longer  does  the  synoptic-scale  environment  indepen¬ 
dently  modulate  the  convection  but  is  itself  modified.  The  degree  of  modification  has  been 
determined  through  detailed  analysis  of  satellite  imagery  and  atmospheric  soundings. 

1.1  Observational  Evidence 

Close  examination  of  satellite  loops  and  upper-air  data  indicates  that  mass  injected 
into  the  upper-troposphere  and  lower-stratosphere  by  convection  causes  systematic  varia¬ 
tions  in  the  flow.  By  computing  the  displacement  of  cirrus  cloud  elements  atop  a  tornado- 
producing  thunderstorm,  Ninomiya  (1971a,  b)  showed  that  the  upper- tropospheric  flow 
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evolves  into  a  difluent  jet  with  maximum  winds  on  the  poleward  side  of  the  convection. 
Maddox  et  al.  (1981)  and  Wetzel  et  al.  (1983)  verified  that  the  jet  streak  is  produced  by 
the  MCS  and  not  by  baroclinicity  in  the  synoptic-scale  flow.  Using  streamline  analysis, 
Fankhauser  (1974),  Johnson  et  al.  (1990),  and  Biggerstaff  and  Houze  (1991)  found  that 
the  outflow  of  mature  squall  lines  is  also  strongly  difluent.  Even  more  impressive  is  the 
intense  divergence  observed  in  the  outflow  of  mature  MCCs  (Pritsch  and  Maddox,  1981a, 
hereafter  FM;  Maddox  et  al.,  1981;  Wetzel  et  al.,  1983;  Cotton  et  ah,  1989;  and  Tripoli  and 
Cotton,  1989a,  b).  FM  suggest  that  the  observed  divergence  “strongly  resembles  the  cir¬ 
culation  which  would  result  if  a  long-lived,  meso-a  scale,  high-pressure  perturbation  were 
superimposed  upon  a  mean  westerly  flow.”  In  fact,  earlier  studies  had  revealed  a  positive 
geopotential  anomaly  atop  the  convection  (Fankhauser,  1974).  Further  objective  analysis 
conducted  by  FM  show  that  a  large  anticyclonic  lens  develops  in  the  upper-tropospheric 
flow. 

Examining  the  upper  air  temperature  data,  FM  discovered  that  the  anticyclone  is 
baroclinically  driven  by  a  pool  of  warm  air  beneath  the  anticyclone  and  a  stronger  cold 
pool  above.  The  temperature  perturbations  were  collocated  with  a  layer  of  moisture  above 
600  mb  indicating  the  presence  of  deep  penetrative  convection.  Similar  thermal  features 
to  those  identified  by  FM  were  also  observed  by  Bosart  and  Nielsen  (1993,  hereafter  BN) 
in  their  analysis  of  soundings  from  a  Gulf  Coast  convective  cluster.  They  found  that  the 
lapse  rate  within  the  outflow  layer  is  nearly  adiabar  .c  and  that  the  wet  bulb  potential 
temperature  is  identical  to  that  at  the  surface.  Furthermore,  reduced  wind  speeds  in  this 
layer  indicate  the  presence  of  low  momentum  air  advected  from  lower  levels.  BN  suggest 
that  the  sounding  penetrated  undilute  air  from  an  updraft  that  had  ascended  from  the 
planetary  boundary  layer. 

1,2  Numerical  Evidence 

Numerical  simulations  of  mesoscale  convection  consistently  produce  the  kinematic 
and  thermodynamic  features  discussed  above  (Fritsch  and  Maddox,  1981b;  Fritsch  and 
Brown,  1982;  Tripoli  and  Cotton,  1989a,  b;  Zhang  and  Gao,  1989;  Schubert  et  al.,  1989; 
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Chen,  1990;  Hertenstein  and  Schubert,  1991;  Zhang,  1992;  Davis  and  Weisman,  1994;  01s- 
son,  1994;  Jiang  and  Raymond,  1995).  Results  from  these  simulations  indicate  that  the 
positive  pressure  anomaly  is  generated  hydrostatically  by  a  pool  of  cold  air  which  forms 
atop  the  convection  (Fritsch  and  Maddox,  1981b),  Four  possible  sources  of  this  cold  air 
„ave  been  identified:  1)  adiabatic  cooling  of  the  overshooting  updraft,  2)  adiabatic  cool¬ 
ing  of  ascending  air  produced  by  mesoscale  heating,  3)  the  downward  heat  flux  produced 
by  cloud-top  entrainment,  and  4)  cloud-top  radiative  cooling.  The  resulting  imbalance 
produced  by  the  pressure  anomaly  leads  to  geostrophic  adjustment  and  the  development 
of  the  observed  mesoscale  anticyclone.  Since  the  duration  and  scope  of  the  adjustment  are 
large  in  relation  to  the  inertial  time  scale  and  the  Rossby  radius  of  deformation,  respec¬ 
tively,  Cotton  et  al.  (1989,  hereafter  CLMT)  has  shown  that  MCSs  represent  inertially 
stable  systems  that  are  in  approximate  geostrophic  balance.  This  realization  has  led  to 
the  development  of  balancr  models  to  both  forecast  the  balanced  flow  (Schubert  et  al, 
1989;  Raymond  and  Jiang,  1990;  Hertenstein  and  Schubert,  1991;  Jiang  and  Raymond, 
1995)  and  to  diagnose  the  balanced  flov/  from  primitive  equation  simulations  (Davis  and 
Weisman,  1994;  Olsson,  1994). 

Many  numerical  studies  of  MCCs  (Fritsch  and  Chappell,  1980;  Fritsch  and  Maddox, 
1981b;  and  Fritsch  and  Brown,  1982;  Zang  and  Fritsch,  1986)  and  squall  lines  (Schu¬ 
bert  et  al.,  1989;  Hertenstein  and  Schubert,  1991)  have  been  developed  using  prescribed 
heating  profiles  to  simulate  the  mesoscale  convection;  however,  heating  is  not  the  only 
method  available  to  parameterize  diabatic  effects.  Convective  motions  •: onerated  by  the 
heating  lead  to  the  diabatic  transport  of  mass.  As  a  result,  diabatic  effects  may  also  be 
parameterized  using  the  transfer  of  mass  across  isentropic  surfaces.  One  of  the  first  to 
use  this  method  was  Raymond  (1983)  in  his  study  of  wave-CISK.  Shutts  has  also  used  it 
extensively  in  the  study  of  deep  penetrative  convection  (Shutts,  1987;  Shutts  et  al.,  1988; 
Shutts,  1994;  Shutts  and  Gray,  1994). 
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1.3  Overview 

The  objective  of  this  study  is  to  show,  using  a  balanced  model,  that  the  observed 
anticyclone  and  accompanying  temperature  perturbations  can  be  partially  explained  by 
the  almost  instantaneous  injection  of  undilute  lower-tropospheric  air  into  the  outflow  of 
mesoscale  convection  followed  by  gradient  adjustment.  In  addition,  the  moderating  effects 
of  cloud-top  cooling  will  be  considered.  The  model  is  developed  using  the  axisymmetric, 
quasi-static  primitive  equations  for  inviscid,  adiabatic  flow  on  an  /-plane.  Previous  efforts 
to  study  the  adjustment  problem  analytically  have  concentrated  primarily  on  the  linearized 
Boussinesq  (Gill  et  al.,  1979;  Shutts,  1987;  McDonald,  1990;  McDonald,  1992)  or  anelastic 
systems  (Shutts,  1994;  Shutts  and  Gray,  1994).  This  work  is  more  general  in  that  the  fluid 
is  continuously  stratified,  and  nonlinear  effects  axe  maintained  through  the  application  of 
gradient  balance  and  the  inclusion  of  advection.  In  order  to  obtain  an  invertibility  principle 
similar  to  the  classic  Rossby  (1938)  adjustment  problem,  isentropic  and  potential  radius 
coordinates  (Schubert  and  Hack,  1983)  are  used  to  eliminate  the  advective  processes, 
making  them  implicit  in  the  coordinate  transformation.  Before  discussing  the  model 
in  more  detail  (Chapter  3),  we  will  review  the  literature  relevant  to  the  development 
and  interpretation  of  the  model  and  its  results  (Chapter  2).  A  control  simulation  is 
then  developed  and  compared  ro  the  observations  of  FM  and  BN  followed  by  a  series  of 
sensitivity  studies  designed  to  examine  the  variability  of  the  mesohigh  to  changes  in  the 
forcing,  static  stability,  latitude,  and  cloud-top  cooling  (Chapter  4).  In  the  conclusion, 
we  will  consider  the  systematic  impact  that  these  mesoscale  anticyclonic  perturbations 
have  on  our  large-scale  NWP  model  forecasts  and  speculate  on  future  research  needs 
(Chapter  5). 


Chapter  2 


REVIEW 


2.1  Introduction 

In  the  introduction,  the  impact  of  sustained  deep  penetrative  convection  within  MCSs 
on  the  kinematic  and  thermodynamic  structure  of  the  upper-troposphere  and  lower- 
stratosphere  was  discussed.  Specifically,  an  anticyclonic  circulation  forms  with  cold  and 
warm  temperature  anomalies  above  and  below,  respectively.  Several  mechanisms  have 
been  suggested  as  the  cause  of  this  structure.  This  chapter  will  provide  additional  detail 
to  the  observational  and  numerical  studies  discussed  in  Chapter  1  (Sections  2.2  and  2.3) 
followed  by  a  summary  of  the  suggested  formation  mechanisms  (Section  2.4).  In  addition, 
material  relevant  to  the  gradient  adjustment  problem  (Section  2.5),  and  the  modifying 
effects  of  entrainment  and  radiation  (Section  2.6)  will  also  be  reviewed. 

2.2  Examining  the  Wind  Anomalies 

FM  hoped  to  extract  the  anticyclonic  circulation  from  the  observed  wind  field  by  sub' 
tracting  a  numerically  fox'ecasted  flow  excluding  the  effects  of  deep  convection.  Since  the 
LFM’s^  convective  parameterization  is  an  adjustment  scheme  (not  including  the  vertical 
fluxes  of  mass,  moi:  rure,  and  momentum  by  the  convection),  the  influence  of  mesoscale 
convection  on  the  forecast  wind  fields  should  be  absent.  Thus,  by  taking  the  vector  dif¬ 
ference  of  the  200  mb  observed  and  12  hour  LFM  forecast  winds  for  the  same  time  should 
yield  the  anomalous  upper  tropospheric  flow  caused  by  the  convection.  Notice  that  this 


^The  Limited-area  Fine-mesh  Model(LFM)  is  a  finite  difference,  primitive  equation  model  developed 
by  the  National  Meteorological  Center  (NMC).  It  is  no  longer  used  by  NMC  to  produce  national  forecasts. 
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method  also  provides  a  measur  ;  of  model  error  caused  by  v,ue  weaknesses  in  the  convective 
parameterization  scheme.  FM  applied  this  technique  to  two  large  MCCs  that  formed  over 
the  Midwest  United  States  on  7  May  1978  (Figs.  2. la-2. 3a)  and  27  June  1979  (Figs.  2.1b- 
2.3b). 

FM’s  objective  analysis  reveals  a  region  of  strongly  divergent,  anticyclonic  outflow 
with  an  average  magnitude  of  approximately  15-20  ms”^  (Fig-  2.3a  and  b).  Similar  anti¬ 
cyclonic  circulations  of  5-10  ms”^  have  been  observed  in  the  outflow  of  tropical  cloud  clus¬ 
ters  (Leary,  1979).  In  both  Ccises  considered  by  FM,  the  jet  was  intensified  on  the  poleward 
side  of  the  convection  and  weakened  to  the  south  with  the  strongest  wind  anomalies  located 
downwind  of  the  convection.  Taking  the  vector  difference  of  mesoscale  model  simulations 
including  and  excluding  the  effects  of  convection,  Maddox  et  al.  (1981)  demonstrated  that 
the  jet  streak  is  produced  by  the  mesoscale  convection  and  not  by  the  large-scale  baroclin- 
icity  of  the  polar  jet.  According  to  Wetzel  et  al.  (1983),  the  shallow  depth  of  the  jet  streak 
also  distinguishes  it  from  the  much  deeper  polar  jet.  FM  further  analyzed  the  7  May  1978 
case  using  a  scale  separation  technique  developed  by  Maddox  (1980b). 

In  this  method,  a  low-pass  filter  first  smoothes  the  data  fields  followed  by  a  1500  km 
bandpass  filter  which  removes  the  mesoscale  modes.  The  recovered  fields  represent  the 
mesoscale  and  macroscale  components  of  the  flow.  Fig.  2.4  beautifully  illustrates  the 
anticyclonic  outflow  with  maximum  speeds  of  20  ms“^  to  the  northwest  and  northeast  of 
the  convection.  Later  observational  studies  have  continued  to  increase  our  understanding 
of  these  features. 

Maddox  (1983)  and  CLMT  used  a  composite  analysis  technique  to  construct  profiles 
of  mean  divergence,  vertical  velocity  a;,  and  relative  vorticity  for  various  stages  in  the 
evolution  of  MCCs.  The  profiles  obtained  by  CLMT  are  shown  in  Fig.  2.5.  During  the 
life  cycle  of  the  storm,  the  level  of  maximum  vertical  velocity  climbs  from  650  to  350  mb 
and  the  magnitude  increases  by  about  25%.  The  increased  ascent  of  mass  impinging 
on  the  tropopause  intensifies  divergence  (and  ultimately  the  anticyclonic  vorticity)  by 
approximately  2/3.  Notice  the  strong  surface  convergence  and  surface  cyclonic  flow.  Most 
of  the  initial  convergence  occurs  in  a  layer  below  700  mb  with  a  maximum  at  the  surface 
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Figure  2.1:  200  mb  observed  wind  field  at  1200  UTC  on  (a)  7  May  1978  and  (b) 
27  June  1979.  The  shaded  area  represents  the  cirrus  cloud  tops  with  temperature 
<  -52°C.  A  full  wind  barb  is  5  ms~^  while  a  flag  is  25  ms~^  (From  Fritsch  and  Maddox, 
1981a). 


Figure  2.2:  Same  as  Fig.  2.1  except  for  the  200  mb  12  hour  LFM  forecast  wind  field  (From 
Fritsch  and  Maddox,  1981a). 


Figure  2.3:  Same  as  Fig.  2.1  except  for  the  vector  difference  of  the  200  mb  observed  and 
12  hour  LFM  forecast  wind  fields  (From  Fritsch  and  Maddox,  1981a). 
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f  igure  2.4:  ^.lesoliigh  perturbation  isotachs  (dashed,  5  ms“^  contour  interval)  and  stream- 
iines  (solid)  obtained  using  a  1500  km  bandpass  filter  on  the  1200  UTC  200  mb  observed 
wind  field  from  the  7  May  1978  case.  The  shaded  area  represents  the  cirrus  cloud  tops 
with  temperature  <  — 52°C  (From  Fritsch  and  Maddox,  1981a). 

that  is  nearly  twice  as  strong  as  the  pre-storm  environment.  As  the  system  matures,  the 
convergence  moves  to  the  mid-levels  of  the  storm.  The  low  and  mid-level  structure  within 
MCSs  may  differ  significantly  between  system  types  (i.e.,  squall  lines  vs.  MCCs);  however, 
the  kinematic  and  thermodynamic  structure  atop  the  updraft  is  often  similar  (Zhang  and 
Gao,  1989;  Chen,  1990;  Fritsch  and  Maddox,  1981b;  Zhang,  1992;  Davis  and  Weisman, 
1994;  Olsson,  1994).  It  appears  that  quasi-balanced  anticyclonic  outflow  is  a  common 
feature  of  long-lived  MCSs  regardless  of  mid  and  low-level  structure.  If  thermal  wind 
balance  is  to  be  maintained,  the  mesohigh  must  be  accompanied  by  unique  temperature 
anomalies. 

2.3  Examining  the  Temperature  Anomalies 

Examining  the  upper  air  data  for  the  7  May  1978  MCC,  FM  detected  anomalously 
warm  and  cold  pools  of  air  in  the  upper  troposphere.  Fig.  2.6  shows  the  150  and  300  mb 
temperature  analysis.  A  6-8°C  intense  cold  pool  is  present  at  150  mb  with  a  weaker  1-2°C 
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Figure  2.5:  MCC  (a)  horizontal  divergence,  (b)  vertical  velocity  (w),  and  (c)  vertical 
component  of  relative  vorticity  profiles  for  the  MCC-12  h,  initial,  mature,  and  dissipation 
stages  (From  Cotton,  et  al.,  1989). 
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warm  pool  below  at  300  mb.  In  addition,  FM  point  out  that  the  soundings  show  a  layer  of 
moisture  beginning  at  about  600  mb  indicating  the  presence  of  deep  convection.  Similar 
thermal  features  to  those  identified  by  FM  were  also  observed  by  BN  in  their  analysis  of 
soundings  from  a  Gulf  Coast  convective  cluster  that  formed  on  26  April  1991. 

The  sounding  of  interest  to  BN,  launched  from  Lake  Charles,  Louisiana  (LCH),  pen¬ 
etrated  the  western  edge  of  the  cluster’s  anvil.  Fig.  2.7  shows  the  soundings  for  LCH 
and  surrounding  stations.  The  conditions  at  Jackson,  Mississippi  (JAN)  and  Slidell, 
Louisiana  (SIL)  are  assumed  to  represent  the  undisturbed  environment  and  adjusting 
environment,  respectively.  Close  inspection  of  the  LCH  sounding  reveals  the  develop¬ 
ment  of  a  nearly  adiabatic  layer  between  185  and  135  mb.  Danielsen  (1982)  has  observed 
similar  adiabatic  layers  within  the  outflow  of  tropical  mesoscale  convection.  Notice  that 
the  wet  bulb  potential  temperature  for  this  layer  is  identical  to  that  at  the  surface.  Re¬ 
duced  wind  speeds  in  this  layer  (as  compared  to  those  at  SIL)  indicate  the  presence  of 
low  momentum  air  advected  from  lower  levels.  BN  suggest  that  the  sounding  penetrated 
undilute  air  from  an  updraft  that  had  ascended  from  the  planetary  boundary  layer.  As 
the  moist,  well-mixed,  potentially  warmer  air  is  injected  into  the  lower-stratosphere  it 
produces  the  adiabatic  layer  with  an  accompanying  stable  layer  between  280  and  250  mb. 
Notice  also  that  the  tropopause  ascends  to  the  top  of  the  adiabatic  layer.  As  shown  in 
the  SIL  sounding,  mixing  causes  the  depth  of  the  adiabatic  and  mixed  layers  to  increase 
downward  while  weakening  their  intensity.  The  next  section  will  summarize  the  proposed 
formation  mechanisms  for  the  kinematic  and  thermodynamics  structures  discussed  above. 

2.4  Explaining  the  Anomalies 

As  the  convective  updraft  overshoots  the  tropopause  and  enters  the  stratosphere,  it 
becomes  significantly  colder  than  its  environment.  According  to  FM,  if  this  source  of  cold 
air  can  be  maintained  by  a  continuous  updraft,  then  a  quasi-permanent  dome  of  cold  air 
develops  with  a  hydrostatically  induced  high  pressure  at  its  base.  The  cold  pool  may 
be  further  enhanced  by  the  net  downward  transport  of  heat  produced  by  the  turbulent 
entrainment  of  warm  stratospheric  air  into  the  outflow  (Shutts  and  Gray,  1994).  This 


Figure  2.6:  Temperature  analysis  at  (a)  150  and  (b)  300  mb  for  1200  UTC  on  7  May  1978. 
Shading  indicates  (a)  — 52°C  or  colder  cloud  tops,  and  (b)  temperatures  warmer  than 
— 34°C  (From  Fritsch  and  Maddox,  1981a). 
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Figure  2.7:  1200  UTC,  26  April  1991  soundings  from  Lake  Charles,  LA  (LCH);  Slidell, 
LA  (SIL);  and  Jackson,  MS  (JAN)  (From  Bosart  and  Nielsen,  1993). 
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effect  will  be  considered  in  greater  detail  later  in  this  chapter  (Section  2.6).  Using  a 
mesoscale  model,  Fritsch  and  Brown  (1982)  simulated  the  cold  pool  by  adding  a  shallow 
Ic /er  of  cooling  at  the  top  of  their  convective  heating  profile.  A  region  of  high  pressure 
was  produced;  however,  it  also  developed  (although  weaker)  when  the  cooling  wais  absent. 
Another  mechanism  was  implicitly  present  in  their  model. 

Ninomiya  (1971a)  and  Fankhauser  (1974)  theorized  that  mesoscale  heating  in  the 
upper  levels  of  the  convection  produces  vertical  motions  and  adiabatic  cooling.  Fritsch 
and  Brown  (1982)  found  that  this  source  of  cold  air  could  also  produce  a  region  of  high 
pressure  albeit  weaker  than  that  produced  by  the  updraft,  Wetzel  et  al.  (1983)  propose 
that  the  heating  increases  the  thickness  of  the  layer,  producing  the  mesohigh.  Their 
observations  show  an  increase  of  60  m  in  the  layer  thickness.  As  the  ascending  mass  from 
both  the  updraft  and  heating  encounters  the  tropopause,  it  loses  its  buoyancy  and  spreads 
laterally. 

If  the  injection  of  mass  in  the  upper  troposphere  and  low  r  stratosphere  continues  for  a 
period  longer  than  the  inertial  time  scale  (/”^),  the  flow  will  geostrophically  adjust  toward 
gradient  balance  forming  an  anticyclonic  lens.  The  amount  of  adjustment  will  depend 
on  the  horizontal  scale  of  the  outflow  in  relation  to  the  Rossby  radius  of  deformation 
(A  =  Cn//).  Thermal  wind  balance  requires  that  this  circulation  be  accompanied  by  a 
cold  (warm)  temperature  anomaly  above  (below).  In  Section  2.5,  we  will  discuss  the 
gradient  adjustment  process  in  more  detail.  Recall  BN’s  observation  that  the  lens  initially 
contains  less  momentum  than  that  of  the  tropospheric  environment. 

A  lens  of  low  momentum  air  at  the  tropopause  may  act  as  an  obstacle  to  the  mean 
flow  (FM;  Wetzel  et  al.,  1983).  Bernoulli’s  principle  states  that  a  decrease  in  speed  as  the 
mean  flow  encounters  the  obstacle  will  result  in  an  increase  in  pressure.  This  is  the  least 
likely  of  the  suggested  mechanisms  since  the  observed  high  pressure  is  located  above  the 
convection  and  not  upstream.  A  final  mechanism  offered  by  Fritsch  and  Brown  (1982), 
indicates  the  importance  of  radiation. 

Cirrus  clouds  undergo  radiational  cooling  to  space  along  their  tops  and  radiational 
warming  from  the  surface  along  their  bases.  Such  a  gradient  in  the  radiation  profile  may 
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strengthen  the  dynamically-forced  cold  and  warm  anomalies  mentioned  previously  (FM; 
Fritsch  and  Brown,  1982).  This  effect  will  be  discussed  in  more  detail  in  Section  2.6,  In 
the  next  section,  we  will  consider  the  impact  of  gradient  adjustment  on  the  balanced  state 
of  the  system. 

2.5  Understanding  Gradient  Adjustment 

Most  of  the  early  investigations  of  atmospheric  convection  and  gradient  adjustment 
concentrated  on  the  formation  of  cyclonic  vortices  by  the  compensating  convergence  of 
mass  and  angular  momentum  following  the  accent  of  convective  thermals  (Inman,  1966; 
McIntyre,  1967;  Sasaki  and  Friday,  1967;  Wilkins  et  ah,  1969;  Friday,  1967,  1969;  McDon¬ 
ald,  1990,  1992).  Little  attention  wats  given  to  the  circulations  that  form  as  the  convection 
comes  to  rest  at  the  level  of  neutral  buoyancy.  The  earliest  examination  of  mass  injections 
and  gradient  adjustment  came  not  from  atmospheric  studies  but  from  studies  of  ocean 
dynamics. 

Observations  have  revealed  the  presence  of  anticyclonic  eddies  within  the  oceans  at  a 
variety  of  scales  (Armi  and  Zenk,  1984;  Manley  and  Hunkins,  1985;  Dugan  et  ah,  1982). 
Oceanographers  have  speculated  that  some  of  these  vortices  are  caused  by  deep  water 
formation  (Killworth,  1976,  1983;  Gascard,  1978)  and  hydrothermal  venting  (Baker  et 
al.  1989;  Speer,  1989).  Both  numerical  (Crepon  et  ah,  1989;  Barnier,  et  al.,  1989;  Mc¬ 
Donald,  1992)  and  laboratory  studies  (Helfrich,  1994)  have  been  conducted  to  understand 
these  phenomena.  Most  notably,  Gill  et  al.  (1979)  and  Gill  (1981)  used  numerical  and 
laboratory  experiments  to  understaiiu  the  adjusted  state  that  develops  following  the  for¬ 
mation  of  deep  water  in  the  Greenland  Sea. 

Using  a  homogeneous,  two-layer,  shallow  water  model.  Gill  et  al.  (1979)  examined  the 
flow  that  develops  following  the  gradual  transfer  of  fluid  from  the  upper  to  the  lower  layer. 
Model  calculations  clearly  show  the  interface  doming,  or  geopotential  change,  between  the 
two  fluids  with  a  corresponding  geostrophically  balanced,  anticyclonic  circulation.  Similar 
features  to  those  obtained  by  the  model  were  observed  in  laboratory  experiments  with  the 
exception  of  cyclonic  vorticity  near  the  dome  center  associated  with  momentum  transfer 
from  the  cyclonically  rotating  mass  sink. 
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Expanding  on  the  two-layer  model,  Gill  (1981)  studied  the  instantaneous  injection  of 
a  homogeneous  layer  of  fluid  into  a  uniformly  stratifled,  rotating  environment.  Assuming 
that  the  injection  is  instantaneous  reduces  the  solutioi.  to  an  initial  value  problem.  Also, 
by  making  the  injection  homogeneous  and  transforming  to  semi-geostrophic  coordinates, 
the  solution  is  isolated  to  the  exterior  fluid  in  the  tovm  of  a  linear  Laplace  equation.  Gill 
found  that  the  solution  was  not  dependent  on  the  transient  state  but  only  on  the  ini¬ 
tial  potential  vorticity  distributions.  Fig.  2.8  shows  the  model  derived  temperature  and 
velocity  distribution.  Notice  the  lenticular  structure  in  the  temperature  field  with  corre¬ 
sponding  baroclinic  anticyclone.  The  aspect  ratio  of  the  lens  scales  as  A/F?  =  N/f^  where 
A  and  H  are  the  Rossby  deformation  radius  and  depth,  respectively,  and  N  and  /  are  the 
buoyancy  and  inertial  frequencies,  respectively.  Numerical  experiments  including  viscous 
effects  indicate  that  horizontal  diffusion  spreads  the  temperature  field  laterally,  causing 
the  lens  to  flatten  and  the  circulation  to  weaken.  Laboratory  experiments  conducted  by 
Hedstrom  and  Armi  (1988)  have  verified  the  viscous  dependent  aspect  ratio.  In  numerical 
experiments  similar  to  Gill’s,  Shutts  (1987),  Shutts  et  al.  (1988),  Shutts  (1994),  and  Shutts 
and  Gray  (1994)  have  studied  the  impact  of  mass  intrusions  in  the  atmosphere. 

Using  a  linearized,  rectilinear,  Lagrangian,  parcel  model,  Shutts  (1987)  examined 
the  thermodynamic  and  kinematic  structure  of  a  Boussinesq  atmosphere  following  the 
rearrangement  of  the  parcels.  The  balanced  state  has  a  lens  structure  at  the  level  of 
neutral  buoyancy  similar  to  that  obtained  by  Gill  (1981),  with  a  front  along  the  axis  of 
tho  sink  formed  by  the  compensating  convergence.  The  aspect  ratio  of  the  lens  is  consistent 
with  Gill’s  formulation.  Shutts  et  al.  (1988)  extended  the  model  to  the  axisymmetric  case 
using  semi-geostrophic  coordinates. 

Fig.  2.9  shows  the  results  of  this  model  for  a  mesoscale  convective  simulation.  Notice 
that  the  vertical  surface  front  does  not  form  in  this  case  due  to  the  conservation  of  angular 
momentum;  instead,  a  sloping  front  or  ''eye  wall”  forms  as  mass  is  drawn  down  into  the 
sink.  In  addition  to  the  balanced  dynamics,  Shutts  (1994)  and  Shutts  and  Gray  (1994) 
also  considered  the  transient  aspects  of  the  geostrophic  adjustment  problem. 

Shutts  (1994)  showed  that  the  effect  of  rotation  is  to  trap  a  portion  of  the  energy 
in  the  balanced  flow  that  would  have  otherwise  been  radiated  away  as  gravity-inertia 
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Figure  2.8:  Contours  of  temperature  (right)  and  velocity  (left)  from  an  analytic  model  of 
a  homogeneous  intrusion  of  fluid  into  a  rotating,  continuously  stratified  environment.  The 
contour  interval  on  both  sides  is  0.1  nondimensional  units.  The  velocities  scale  as  NH 
(From  Gill,  1981). 

waves.  The  amount  of  energy  trapped  is  dependent  on  the  scale  of  the  disturbance  in 
elation  to  A.  Using  a  rectilinear,  anelastic,  moist  convective  plume  model,  Shutts  and 
Gray  (1994)  showed  that  30-40%  of  the  initial  APE  (APEq)  is  partitioned  to  the  final 
balanced  kinetic  energy  (KE).  This  is  approximately  nn  order  of  magnitude  greater  than 
that  predicted  by  Schubert  et  al.’s  (1980)  linear  analysis  (Shutts  and  Gray,  1994).  Hack 
and  Schubert  (1986)  showed  that  the  amount  of  energy  partitioned  to  the  balanced  flow 
increases  as  the  inertial  stability  of  the  vortex  increases.  Partitioning  the  flow  from  an 
MCC  simulation  into  balanced  and  unbalanced  components,  Olsson  (1994)  found  that 
during  the  storm’s  life,  as  much  as  84%  of  the  total  KE  may  be  partitioned  into  the 
balanced  KE. 

CLMT  explain  that  the  scale  of  mesoscale  disturbances  is  determined  not  by  the  linear 
form  of  the  Rossby  radius  of  deformation  A  =  c„//  but  by  the  general  form  (Schubert  et 
aJ.,1980) 


Figure  2.9:  Results  from  an  axisymmetric,  Lagrangian  parcel  model.  Fields  shown  are  the 
(a)  initial  element  configuration  with  shaded  area  being  the  positively  buoyant  elements 
(2. IK  increase  in  6  with  each  element),  (b)  final  element  configuration  with  lens  and  “eye” 
wall  structure,  (c)  parcel  displacement  field,  and  (d)  the  balanced  tangential  velocity  (1 
ms“^  contour  interval).  Initially,  elements  in  the  same  horizontal  plane  have  the  same 
potential  temperature  while  elements  in  the  same  vertical  column  have  the  same  angular 
momentum  (From  Shutts  et  al.,  1988). 
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A'  = 


+  (2;^  +  /) 


(2.1) 


where  Cn  is  the  phase  speed  of  a  gravity-inertia  wave,  C  is  the  vertical  component  of  relative 
vorticity,  /  is  the  Coriolis  parameter,  Pic  is  the  radius  of  curvature,  and  v  is  the  tangential 
wind  component  at  iJc-  If  the  convection  is  distributed  over  a  scale  smaller  than  A',  most 
of  the  available  potential  energy  (APE)  generated  by  the  convection  will  be  radiated  away 
as  gravity-inertia  waves.  On  the  other  hand,  if  the  convection  is  distributed  over  a  scale 
larger  than  A',  most  of  the  APE  will  be  imparted  to  the  balanced  flow.  Notice  that  A'  is 
dependent  on  the  relative  vorticity  of  the  system  such  that  an  increase  of  relative  vorticity 
leads  to  a  decrease  of  A^  Thus,  the  scale  of  the  system  changes  as  it  develops  maJ^ing  the 
terms  large  and  small  relative.  Emphasizing  this  point,  CLMT  state  that  a  “disturbance 
having  a  radius  of  700  km  in  the  tropics  may  be  dynamically  small,  while  a  similar-sized 
system  in  midlatitudes  would  be  classified  as  dynamically  large.” 

Defining  the  radius  of  curvature  Rc  as  the  approximate  radius  of  the  area  covered 
by  the  isotherm  .  CLMT  fouiic  :hat  for  their  study  Rc  ~  322  km.  Furthermore, 

using  a  value  of  Cn  ^  30ms”  they  found  that  A'  «  300  km  with  an  accuracy  of  ±4%. 
CLMT  inferred  from  this  that  MCCs  are  actually  inertially  stable  systems.  This  result 
will  be  important  when  defining  the  horizontal  magnitude  of  mass  perturbations  and 
evaluating  the  final  balanced  state  of  the  model.  In  the  next  section,  the  effects  of  cloud- 
top  entrainment  and  radiative  cooling  on  the  magnitude  of  the  cold  pool  will  be  discussed. 


2.6  Revealing  the  Effects  of  Entrainment  and  Radiation 


2.6.1  Cooling  by  Entrainment 

When  shallow  convection  encounters  an  inversion,  it  entrains  warm  air  from  above  the 
inversion.  Ball  (1960)  showed  that  this  entrainment  results  in  a  net  downward  transport 
of  heat,  causing  the  height  of  the  inversion  to  increase  while  the  layer  below  remains  well 
mixed  and  nearly  moist  adiabatic.  This  process  is  essential  to  the  dissipation  of  turbulent 
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kinetic  energy  and  the  production  of  potential  energy.  In  his  studies  of  stratosphere- 
troposphere  exchange,  Danielsen  (1968,  1982,  1993)  has  proposed  a  similar  m<  anism 
(shown  in  Fig.  2.10)  for  deep  penetrative  convection  encountering  the  tropopause. 


Figure  2.10:  Theoretical  mechanism  for  stratospheric  anvil  formation  and  its  associated 
temperature  distribution  (From  Danielsen,  1982). 

A  ’  ieep  convection  impinges  on  the  tropopause,  updrafts  overshoot  the  level  of  neutral 
buoyancy  and  enter  the  stratosphere.  On  average,  vertical  velocities  are  about  20-60  ms“^ 
with  penetrations  of  1-3  km  (Saunders,  1962).  Increased  static  stability  within  the  strato¬ 
sphere  causes  the  updraft  to  lose  buoyancy  and  decelerate.  Like  a  giant  fountain,  the 
decelerating  air  eventually  collapses  down  onto  the  updraft. 

Turbulent  eddies,  associated  with  shear  instabilities,  develop  at  the  boundary  of  the 
overshooting  updraft  (Klaassen  and  Clark,  1985;  Clark  et  al.,  1988;  Grabowski,  1989; 
Grabowski  and  Clark,  1991,  1993)  and  entrain  warm  stratospheric  air.  If  the  level  of 
neutral  buoyancy  is  the  tropopause,  the  increased  buoyancy  of  the  diluted  air  causes  the 
collapsing  updraft  to  come  to  rest  within  the  stratosphere.  IR  satellite  imagery  (Adler 
and  Fenn,  1979;  Fujita,  1982;  Heymsfield  and  Blackmer,  1988;  Ebert  and  Holland,  1992), 
aircraft  observations  (Adler  et  al.,  1983;  Spinhirne  et  al.,  1983;  Di:  lielsen,  1981,  1993; 
Poulida  et  al.,  1995)  and  numerical  simulations  (Adler  and  Mack,  1986;  Schlesinger,  1990; 
Stenchikov  et  al.,  1995)  support  the  development  of  anvils  within  the  stratosphere. 
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The  entrainment  and  mixing  within  the  stratospheric  outflow  result  in  a  net  down¬ 
ward  transport  of  heat  and  the  production  of  a  well  mixed,  adiabatic  layer  with  sharp 
temperature  gradients  along  its  horizontal  boundaries.  Ackerman  et  al.  (1988)  warns  that 
the  resulting  temperature  profile  may  be  less  than  adiabatic  (or  even  isothermal)  depend¬ 
ing  on  the  structure  and  vigor  of  the  turbulence  within  the  layer.  Notice  that  the  shape  of 
the  LCH  sounding  obtained  by  BN  agrees  well  with  Danielsen’s  theoretical  profile  (com¬ 
pare  Figs.  2.7  and  2.10).  According  to  Danielsen  (1982),  the  mixing  within  the  layer  is 
enhanced  by  the  destabilization  caused  by  differential  radiative  heating  along  the  top  and 
bottom  of  the  anvil. 

.6,2  Cooling  by  Radiation 

Longwave  cooling  along  the  top  of  the  outflow  causes  ^-surfaces  to  “lift” ,  transferring 
mass  to  lower  isentropic  layers  (Wirth,  1995).  These  diabatic  mass  transfer  mechanisms 
must  be  accounted  for  when  considering  the  adjustment  process.  Johnson  et  al.  (1990),  in 
their  analysis  of  squall  lines,  showed  that  diabatic  cooling  by  radiation  leads  to  mesoscale 
descent  along  the  upper  surface  of  the  stratiform  outflow. 

Webster  and  Stephens  (1980)  have  measur  ed  cloud-top  cooling  rates  within  the  trop¬ 
ics  of  approximately  5-15°Cday”^  (~  .5°Chr^^).  Note  that  their  work  concentrates  pri¬ 
marily  on  mid-l  vel  clouds  and  not  on  upper-tropospheric  cirrus.  Simulations  of  tropical 
MCSs  by  Wong  et  al.  (1993  a,  b)  show  cooling  rates  between  2  and  8°Cday~  ' .  Although 
no  current  measurements  exist,  similar  cooling  rates  may  be  expected  atop  mid-latitude 
MCSs. 

Chen  and  Cotton  (1988)  conducted  sensitivity  studies  of  a  High  Plains  MCS  simula¬ 
tion  and  found  2®C  warming  of  the  cloud-top  in  the  absence  of  longwave  radiative  effects. 
Similar  sensitivity  analysis  by  Tao  et  al.  (1993)  for  a  mid-latitude  squall  line  show  cooling 
rates  of  2-3°Cday“‘^. 

Starr  and  Cox’s  (1985a,  b)  emphasize  that  cloud-top  cooling  is  enhanced  during  the 
evening.  In  particular,  the  additional  cooling  could  have  a  significant  impact  on  nocturnal 
systems  such  as  MCCs  (Wetzel  et  ah,  1983).  MCS  simulations  conducted  by  Tripoli  and 
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Cotton  (1989b)  show  an  increase  in  cloud-top  cooling  from  4  to  IG'^Cday  ^  during  the 
evening. 


Chapter  3 


MODEL 


3.1  Introduction 

In  the  previous  chapter,  we  discussed  in  detail  the  complex  kinematic  and  thermo¬ 
dynamic  structure  within  the  lower-stratospheric  outflow  of  an  MCS.  In  addition,  we 
considered  the  moderating  effects  of  rotation,  entrainment,  and  radiation.  As  a  convec¬ 
tive  tower  stretches  through  the  depth  of  the  troposphere,  mass  is  rapidly  redistributed 
producing  an  imbalance  in  the  forces  acting  on  the  fluid.  The  nonlinear  adjustment  process 
that  follows  is  partitioned  into  a  transient  and  balanced  component.  During  the  transient 
adjustment,  excess  mass  and  energy  are  shed  into  the  far-field  through  the  propagation 
of  gravity-inertia  waves,  leaving  behind  a  mass  field  in  hydrostatic  and  gradient  balance. 
Solutions  to  the  nonlinear  transient  adjustment  are  complicated  (Shutts  and  Gray,  1994); 
however,  the  linearized  form  has  been  solved  analytically  (Schubert  et  al.,  1980;  Fulton 
and  Schubert,  1980;  Shutts,  1994).  The  purpose  of  this  chapter  is  to  describe  the  analytic 
development  of  an  axisymmetric,  nonlinear,  balance  model  designed  to  find  the  adjusted 
state  of  an  impulsively  perturbed  system  using  the  conservation  of  potential  vorticity 
(PV). 

Since  the  core  of  this  model’s  formulation  is  based  on  PV  conservation,  it  will  be 
profitable  to  first  review  the  definition  of  PV,  how  convection  generates  PV  anomalies, 
and  how  the  balanced  flow  is  extracted  from  the  PV  field  (Section  3.2).  Given  this  infor¬ 
mation,  we  will  define  the  basic  state  with  appropriate  governing  equations  (Section  3.3) 
and  develop  a  suitable  model  to  extract  the  balanced  state  (Section  3.4).  It  will  also  be 
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beneficial  to  analyze  the  model’s  energy  and  mass  bude::i:s  to  ensure  that  they  are  consis¬ 
tent  (Section  3.5).  The  remainder  of  the  chapter  will  explain  the  numerics  used  in  solving 
the  model  (Section  3.6). 

3.2  Thinking  in  Terms  of  PV 
3,2.1  PV  Definition 

The  foundations  of  PV  thinking  originate  with  the  early  formulations  of  circulation 
and  vorticity  by  Helmholtz  (1858),  Kelvin  (Thompson  1869,  1879),  V.  Bjerknes  (1898), 
and  others.  As  pointed  out  by  Pedlosky  (1979,  §  2.5),  the  vorticity  equation  can  describe 
the  evolution  of  atmospheric  motions  but  places  no  constraints  on  that  evolution;  whereas, 
we  know  that  gravity  (i.e.,  static  stability)  and  planetary  rotation  (i.e.,  inertial  stability) 
place  unique  restrictions  on  fluid  motions  within  the  atmosphere.  In  particular,  thermal 
wind  balance  must  be  maintained. 

As  an  example  of  this  point,  consider  two  axisymmetric,  diabatically  generated  mass 
perturbations  within  a  resting,  statically  stable  atmosphere  as  shown  in  Fig.  3.1a  and 
c.  The  APE  of  the  system  is  increased  by  the  adiabatic  displacement  of  mass  at  the 
layer  of  injection/removal.  As  gravity-inertia  waves  transport  mass  radially,  the  APEq  is 
converted  into  KE.  Ultimately,  the  deflection  of  the  radial  flow  by  the  Coriolis  force  draws 
the  systems  into  geostrophic  balance.  The  resulting  balanced  states  are  shown  in  Fig.  3.1b 
and  d.  Notice  that  a  portion  of  the  APEq  must  remain  in  order  to  maintain  thermal  wind 
balance.  As  stated  by  Shutts  (1994),  “rotation  prevents  the  complete  conversion  of  PE  to 
KE  through  the  conservation  of  angular  momentum.”  Apparently,  a  unique  relationship 
exists  between  changes  in  the  APE  (i.e.,  static  stability)  and  the  absolute  vorticity. 

While  developing  a  new  hydrodynamic  vorticity  equation,  Ertel  (1942)  showed  that, 
for  adiabatic,  inviscid  motions,  the  static  stability  and  absolute  vorticity  are  related  by 
the  following  material  conservation  relation: 


^K«-V0)  =  o, 


(3.1) 
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Figure  3.1:  (a)  Atmosphere  perturbed  by  a  mass  source  “+”  with  (b)  the  resulting  bal¬ 
anced  state,  (c)  Atmosphere  perturbed  by  a  mass  sink  with  (d)  the  resulting  balanced 
state.  The  arrows  indicate  the  direction  of  gravity  wave  propagation  while  the  circle  with 
a  dot  in  the  center  denotes  flow  out  of  the  paper  and  the  circle  with  an  “x”  in  the  center 
denotes  flow  into  the  paper. 


where  a  =  l/p  is  the  specific  volume,  (^a  is  the  three-dimensional  absolute  vorticity,  and  V 
is  the  three-dimensional  gradient  operator.  Amazingly,  no  approximations  are  involved  in 
this  result.  It  is  valid  for  three-dimensional,  nonhydrostatic,  baroclinic  flows  at  any  scale. 
Rossby  (1940)  defined  the  parenthetical  quantity  the  potential  vorticity  (PV) 


P  =  QfCa  •  =  constant.  (3.2) 

Thus,  variations  in  a  parcel’s  static  stability  and  inertial  stability  are  constrained  by  the 
requirement  that  PV  be  materially  conserved.  Again,  as  Shutts  (1994)  explains,  “the  total 
energy  is  minimized  subject  to  the  Lagrangian  conservation  of  potential  temperature  and 
angular  momentum  (Shutts  and  Cullen  1987);  or,  alternatively,  the  conservation  of  PV 
(Vallis  1992).  In  the  non-rotational  case,  energy  minimization  results  in  a  horizontally 
homogeneous  end  state  with  all  the  KE  appearing  in  gravity-wave  motion,”  Equation  (3.2) 
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defines  the  link  between  PV  anomalies,  and  the  kinematic  and  thermodynamic  states  of 
a  system,  but  does  not  address  how  matss  perturbations  generate  PV  anomalies. 

3,2.2  Impermeability  Theorems 

In  a  controversial^  pair  of  papers  by  Haynes  and  McIntyre  (1987,  1990,  hereafter  HM), 
a  set  of  constraints  known  as  the  impermeability  theorems  were  developed  that  set  limits 
on  the  transport  of  PV  and  provide  a  useful  framework  for  understanding  the  development 
of  PV  anomalies.  HM  suggest  that  the  conservation  of  PV  constitutes  a  partial  analogy 
with  the  conservation  of  a  chemical  mixing  ratio.  In  the  absence  of  chemical  reactions, 
diffusion  across  boundaries,  and  sedimentation,  the  mixing  ratio  of  a  chemical  substance 
is  materially  conserved;  however,  unlike  the  chemical  mixing  ratio,  which  is  composed 
of  a  chemical  substance,  PV  is  a  mathematical  construct  with  no  physical  analog  to  the 
chemical  substance.  To  alleviate  this  discrepancy,  HM  define  a  quantity  called  PV  sub- 
stance  (PVS)  recognizing  that,  unlike  the  chemical  substance,  PVS  is  a  signed  quantity. 
HM  demonstrated  mathematically  that  PVS  is  trapped  between  isentropic  surfaces  and 
can  only  be  fluxed  along  the  surfaces.  Based  on  this  result,  they  developed  the  following 
theorems: 

(i)  There  can  be  no  net  transport  of  PVS  across  any  isentropic  surface. 

(ii)  PVS  can  neither  be  created  nor  destroyed,  within  a  layer  bounded  by  two  isentropic 
surfaces.  Corollary:  PVS  can  only  be  created  or  destroyed  where  the  layer  terminates 
along  a  boundary. 

Remarkably,  this  is  true  even  in  the  presence  of  diabatic  and  viscous  processes.  Consider 
how  these  restrictions,  in  conjunction  with  the  mixing  ratio  analogy,  might  be  used  to 
visualize  the  generation  of  PV  anomolies. 


^  Based  on  observational  evidence,  Danielson  (1990)  disputed  Haynes  and  McIntyre’s  (1987)  imperme¬ 
ability  conclusions.  In  a  follow-up  paper,  Haynes  and  McIntyre  (1990)  clarified  their  theory  and  demon¬ 
strated  that  the  differences  were  primarily  seinautic. 
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Fig.  3.2  depicts  the  development  of  two  mass  perturbations  and  their  corresponding 
PV  anomalies  induced  by  the  diabatic  redistribution  of  mass  by  moist  convection.  Let’s 
examine  the  lower  perturbation  first.  As  mass  is  diabatically  removed  from  the  layer,  the 
isentropes  are  drawn  toward  the  perturbation  by  the  adiabatic  descent  of  mass  filling  the 
void  left  by  the  mass  sink.  According  to  theorem  (i),  the  increased  gradient  of  isentropic 
surfaces  concentrates  the  PVS  and  increases  lie  PVS  mixing  ratio  (i.e.,  generating  a  pos¬ 
itive  PV  anomaly).  The  opposite  holds  for  the  upper  perturbation.  Mass  is  diabatically 
injected  into  the  layer  causing  the  isentropes  to  spread  apart.  The  decreased  isentropic 
gradient  dilutes  the  PVS  and  decreases  the  PVS  mixing  ratio  (i.e.,  generating  a  negative 
PV  anomaly).  By  comparing  Figs.  3.1  and  3.2,  we  conclude  that  a  positive(negative)  PV 
anomaly  is  associated  with  cyclonic(anticyclonic)  vorticity.  An  understanding  of  PV  con¬ 
servation  combined  with  the  constraints  imposed  by  the  impermeability  theorems,  provide 
a  qualitative  picture  of  the  formation  and  influence  of  PV  anomalies.  This  information 
alone,  however,  does  not  allow  the  mathematical  extraction  of  the  balanced  flow  from  the 
distribution  of  PV. 

3,2.3  Invertibility  Requirements 

Given  the  perturbed  wind  and  mass  fields,  the  initial  PV  distribution  is  obtainable 
from  (3.2)  while  its  future  distribution  may  be  forecast  by  (3.1);  however,  a  problem 
develops  when  one  attempts  to  extract  the  future  state  of  the  wind  and  mass  fields  from  the 
forecast  distribution  of  PV.  For  a  given  value  of  PV,  (3.2)  indicates  an  infinite  number  of 
kinematic  and  thermodynamic  combinations.  Obviously,  additional  information  is  needed. 

As  pointed  out  by  Hoskins  et  al.  (1985,  hereafter  HMR),  the  following  conditions 
must  be  specified  in  order  to  invert  the  potential  vorticity  field  and  obtain  a  velocity  and 
temperature  field  consistent  with  thermal  wind  balance: 

(i)  A  reference  state  must  be  defined, 

(ii)  Two  balanced  conditions  must  be  specified. 

(iii)  The  inversion  problem  must  be  solved  globally  with  appropriate  boundary  conditions 


given. 
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Figure  3.2:  Generation  of  and  PV  anomalies  by  moist  convection.  Horizontal 
lines  represent  isentropic  surfaces. 

Since  an  anomaly  refers  to  a  deviation  from  some  reference,  condition  (i)  is  necessary 
to  clarify  what  is  meant  by  PV  anomaly.  Conditions  (ii)  and  (iii)  will  allow  us  to  find 
a  unique  balanced  state  for  a  given  distribution  of  PV.  Notice  in  (3.2)  that  the  PV  is 
a  function  of  three  3  variables,  a,  and  V6.  By  specifying  two  balanced  conditions 
relating  the  three  variables,  the  PV  is  reduced  to  a  function  of  one  variable.  The  resulting 
system  of  equations  is  commonly  referred  to  as  the  invertibility  principle.  Since  the  PV 
anomaly  affects  the  flow  both  locally  and  into  the  far-field,  the  invertibility  principle  must 
be  solved  globally  with  appropriate  boundary  conditions  specified.  Thus,  given  an  initial 
distribution  of  PV  and  knowing  that  it  is  materially  conserved,  we  can  solve  for  the  final 
state  of  the  system  consistent  with  the  balance  conditions.  If  a  local  conservation  relation 
can  be  developed,  then  the  solution  reduces  to  an  initial  value  problem.  Each  of  these 
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requirements  will  be  dealt  with  in  following  sections  as  the  model  evolves  beginning  with 
condition  (i)  —  specification  of  the  basic  state. 

3.3  Formulating  the  Problem  of  Moist  Convection  using  PV  Anomalies 

3.3.1  Basic  State  and  Governing  Equations 

For  simplicity,  we  consider  a  stably  stratified  atmosphere  at  rest  relative  to  an  /-plane. 
Motions  within  the  atmosphere  are  assumed  quasi-static,  adiabatic  and  inviscid.  Since 
the  basic  state  is  statically  stable  (i.e.,  6  is  monotonically  increasing  with  height)  and  the 
motions  are  adiabatic,  it  will  be  convenient  to  use  isentropic  coordinates  in  the  vertical. 
However,  instead  of  0,  specific  entropy 

will  be  used,  where  Cp  is  the  heat  capacity  of  dry  air  at  constant  pressure,  Rd  is  the  gas 
constant  of  dry  air,  “o”  denotes  a  constant  reference  value  with  =  To,  and  the  pressure 
p,  density  p  and  absolute  temperature  T  are  related  by  the  ideal  gas  law  p  =  pRdT. 
The  use  of  s  yields  a  more  convenient  form  of  the  hydrostatic  equation  in  which  the 
right  hand  side  is  simply  T  rather  than  the  Exner  function^.  Fig.  3.3a  illustrates  the 
model  basic  state  absolute  temperature  T{s)  profile.  In  order  to  study  the  impact  of 
stratospheric  penetrations  by  convection,  the  tropopause  is  specified  within  the  interior  of 
the  domain  rather  than  using  the  top  boundary.  T  was  chosen  because  of  its  familiarity 
(e.g.,  atmospheric  soundings),  but  it  is  the  density  profile  that  we  will  be  most  interested 
in  since  perturbations  will  be  to  the  mass  field. 

Previous  models  have  focused  on  basic  state  density  profiles  that  were  either  Boussi- 
nesq  (Gill  et  al.,  1979;  Shutts,  1987;  McDonald,  1990;  McDonald,  1992)  or  anelastic 
(Shutts,  1994;  Shutts  and  Gray,  1994).  In  this  model,  we  will  consider  mass  perturba¬ 
tions  superimposed  on  a  continuously  stratified  atmosphere.  Recall  that  the  density  in 


Exner  function,  11  =  Cp(p/po)^- 
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s  s 


Figure  3,3:  Basic  state  (a)  temperature  T{s)  and  (b)  pseudodensity  (t(s)  profiles.  Sstr  is 
the  specific  entropy  difference  over  which  the  tropopause  discontinuity  is  smoothed. 

isentropic  coordinates  is  given  by  the  pseudodensity  cr,  where  pdz  =  —g'^^dp  =  g^^ads 
su that  a  =  —dplds.  The  basic  state  pseudodensity  a{s)  profile  (Fig.  3.3b)  is  derived 
from  T{s)  by  solving  for  the  pressure  and  then  differentiating  with  respect  to  s.  As  an 
additional  simplifying  assumption,  we  will  consider  the  axisymmetric  case.  The  follow¬ 
ing  set  of  primitive  equations,  in  cylindrical  coordinates,  governs  the  motions  within  our 
hypothetical  atmosphere: 


Du  ( .  v\  dM 

(3.4) 

Dv  (  .  v\ 

(3.5) 

II 

(3.6) 

Dg  dim) 

Dt'^^  rdr  ~ 

(3.7) 

where 
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Dt  dt^  dr^ 


(3,8) 


is  the  material  derivative,  u  and  v  are  the  radial  and  tangential  wind  components,  re¬ 
spectively,  and  M  =  CpT  -h  gz  is  the  Montgomery  potential  or  streamfunction.  Poisson’s 
diagnostic  relationship 


r  =  roe"/‘=^  (3.9) 

\pqJ 

may  also  be  added  to  our  equation  set.  Using  the  definition  a  =  —dp/ds,  the  above  system 
of  equations  (3.4)-(3.9)  is  closed  in  the  unknowns  n,  t;,  M,  T,  a,  and  p. 

The  PV  conservation  relation  for  this  system  is  found  by  combining  the  isentropic 
absolute  vorticity  equation 


^  _ 

Dt  rdr 

with  the  continuity  equation  (3.7),  eliminating  the  isentropic  divergence  to  obtain 


(3.10) 


^(7)  =  “’ 

where,  Ca  =  j  +  d{rv)/rdr  is  the  isentropic  absolute  vorticity  and  P  =  Ca/<^  is  the 
PV.  In  later  calculations,  it  will  be  convenient  to  use  the  inverse  of  P  or  the  potential 
pseudodensity  a*  =  f  /P  =  (//Ca)c^  (i-e.,  the  pseudodensity  that  a  parcel  would  have  if 
Ca  =  /)•  With  this  definition,  (3.11)  may  be  rewritten  as 


Da* 

Dt 


=  0. 


(3.12) 


With  the  basic  state  and  governing  equations  defined,  we  need  only  define  how  the  system 
will  be  perturbed. 
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3. 3, 2  Perturbation 

Exactly  what  form  the  perturbation  should  take  and  how  it  should  be  imparted 
on  the  basic  state  presents  some  difRcvilties.  As  discussed  in  Chapter  2,  convection  is 
a  highly  turbulent  process  and  as  such  is  computationally  cumbersome  and  difficult  to 
model;  however,  simplifying  assumptions  to  the  primitive  equations  or  to  the  form  of  the 
perturbation  will  reduce  the  complexity. 

If  the  time  scale  of  the  convection  is  much  less  than  the  time  scale  of  the  adjustment 
(/’■^),  then  the  movement  of  mass  is  assumed  to  occur  instantaneously.  Typically,  the  ver¬ 
tical  motions  within  a  convective  updraft  are  rapid  (10-15  min);  whereas,  the  adjustment 
occurs  over  much  longer  time  scales  (/”^  4  hrs  at  30°N)  (Shutts,  1987;  Helfrich,  1994). 

This  assumption  is  valid  for  most  forms  of  atmospheric  convection.  Even  within  an  MCC, 
where  convection  may  last  for  hours,  the  time  scale  over  which  mass  is  injected  into  the 
upper  troposphere  is  much  less  than  the  adjustment  time  scale  (Maddox,  1983).  By  mak¬ 
ing  this  assumption,  the  complex  turbulent  transport  process  is  excluded.  Observe  that 
the  fluid  could  be  gradually  injected,  preserving  a  state  of  quasi-balance  (Gill  et  al,,  1979). 
Although  useful,  this  method  would  eliminate  our  ability  to  treat  the  problem  as  an  initial 
value  problem.  The  pseudodensity  perturbation  can  be  expressed  as  follows: 


<7'(r-,s,0)  =  aF{r)G{s), 


(3.13) 


where  F{r)  and  G{s)  are  the  horizontal  and  vertical  striicture  functions,  respectively,  and 
a  is  the  amplitude  of  the  perturbation.  The  vertical  and  horizontal  perturbation  profiles 
are  mathematically  given  by 


F(r)  =  { 


0 

l-/[(r-ri)/(r2-ri)] 

1 


r2  <  r  <  oo, 
ri  <  r  <  r2, 
0  <  r  <  ri, 


(3.14) 


and 
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Figure  3.4:  The  (a)  horizontal  F{r)  and  (b)  vertical  G(s)  structure  functions  for  a  hypo¬ 
thetical  MCS. 


where 


'  0 

-eG{l  -  I[(s  -  S4)/(S5  -  S4)]} 

G(s)  =  -  i  -  (1  +  «)-^[(s  -  »3)(S4  -  S3)]} 
2  G/[(s  -  S2){s3  -  S2)] 

0 


S5  <  S  <  ST, 
54  S 
S3  <  S  <  S4, 
S2  <  S  <  S3, 
Si  <  s  <  S2, 
0  <  s  <  Si, 


(3.15) 


/(a;)  =  exp{  — (2x)  Mn(2)  exp[(l  —  2x)/(l  “  x)]}  0<x<l  ,  (3.16) 

is  a  smooth  unit  step  function  with  the  properties  7(0)  =  0,  7(^)  —  and  7(1)  =  1. 
Fig.  3.4  shows  the  structure  functions.  The  constant  G  is  defined  such  that  mass  continuity 
is  maintained  (i.e.,  G{s)ds  =  0) 


G  = 


siil-I) 


(3.17) 


(53  ^  52)7  +  {Sr  —  S^)[l  —  (1  +  e)7]  -f  €(55  —  54)(7  —  1) 
where  7  =  /q^  I{x)dx  ^  0.501407  is  the  mean  of  the  smoothing  function.  Mass  is  removed 
from  a  shallow  surface  layer  0  <  s  <  si  and  injected  into  the  upper-troposphere  or 
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lower-stratosphere  52  <  5  <  54.  Notice  that  this  method  differs  from  that  employed  by 
Gill  (1981)  and  Shutts  (1987)  in  that  the  injected  mass  is  not  homogeneous  but  is  itself 
stratified.  Mass  may  also  be  removed  from  above  the  injection  level  54  <  s  <  S5  simulating 
the  diabatic  effects  of  entrainment  and  radiation.  The  constant  e  determines  the  amount 
of  mass  removed.  Notice  that  if  the  upper  mass  sink  is  included,  a  section  of  the  profile 
within  the  interval  53  <  5  <  S4  becomes  part  of  the  sink.  This  proves  to  be  a  useful 
artifact  of  the  profile  since  the  effects  of  entrainment  and  radiation  do  extend  into  the 
outflow.  Combining  the  initial  values  of  a  and  cr'  yields  the  total  perturbed  mass  field 

a(r,  s,  0)  =  a{s)  +  cr'(r,  s,  0).  (3.18) 

Since  the  atmosphere  is  initially  at  rest,  the  definition  of  cr*  indicates  that 

a^{T,  s,  0)  =  (7(r,  s,  0).  (3.19) 

With  the  initial  distribution  of  cr*  defined,  we  can  develop  a  model  to  extract  the  final 
balanced  state  of  the  system. 

3.4  Extracting  the  Balanced  State  from  the  PV  Distribution 

3.4^1  Gradient  Balance  Condition 

An  instantaneous  perturbation  to  the  mass  field  disturbs  the  horizontal  force  balance 
(/  +  v/t)v  ^  dMfdr  producing  radial  and  tangential  flow  through  (3.4)  and  (3.5),  re¬ 
spectively.  Adjustments  to  the  mass  field  through  continuity  (3.7)  will  gradually  draw  the 
system  toward  a  new  balanced  state. 

HMR’s  provision  (ii)  states  that  two  balanced  conditions  be  specified.  Requiring 
that  motions  be  quasi-static  meets  the  first  condition.  To  maintain  the  highest  level  of 
generality,  we  choose  gradient  balance  as  the  second.  Using  the  definition  for  P,  ^md 
cr,  the  following  set  of  equations  describes  the  balanced  motion: 


(3.20) 
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dM 


(3.21) 


dM 

ds 


(-)"■ 

\PoJ 


(3.22) 


According  to  HMR’s  condition  (iii),  appropriate  boundary  conditions  must  also  be  spec¬ 
ified.  We  shall  assume  that  both  the  upper  and  lower  boundaries  are  isentropic  surfaces. 
Furthermore,  we  shall  require  that  the  upper  boundary  be  an  isobaric  surface  while  the 
lower  boundary  will  be  taken  as  a  constant  height  surface.  Along  the  outer  boundary,  the 
pressure  approaches  a  far-field  value  p{s)  representing  the  mean  mass  distribution.  From 
these  requirements,  our  boundary  conditions  become 


P  =  PT  at  5  =  ST,  (3.23) 

M  =  CpTo  (^)  at  s  =  0,  (3.24) 

p^p  as  r  — ►  oo.  (3.25) 

(3.20) -(3.25)  represents  the  invertibility  principle  for  v(r,  s),  M(r,  s)  and  p(r,  s)  in  terms 
of  a*(r,  s)  and  p(s).  Thus,  given  the  forecast  distribution  of  a*  from  (3.12)  and  p(s), 

(3.20) -(3.22)  are  solved  concurrently  to  find  v,  M  and  p.  Alternatively,  (3.20)-(3.22) 
could  be  combined  into  a  single  second  order  partial  differential  equation  in  terms  of  u, 
M  or  p,  and  a*.  HMR  and  Thorpe  (1985,  1986)  discuss  the  invertibility  principles  for  v 
and  p,  respectively;  however,  we  choose  to  take  the  more  common  approach  that,  in  its 
simplest  form,  inverts  the  Laplacian  of  a  streamfunction  (Raymond,  1992).  In  this  case, 
the  Laplacian  is  a  complicated  second  order  differential  operator  and  the  streamfunction 
is  M.  By  multiplying  (3.21)  by  r^,  differentiating  this  result  with  respect  to  r,  and  using 

(3.20)  and  (3.22)  to  eliminate  v  and  p,  we  obtain  the  invertibility  principle  in  terms  of  M, 
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+  J-  (',3 

r^dr  V  dr  ) 

The  boundary  conditions  (3.24)-(3.25)  in  terms  of  M  become 

at  5  =  sj-,  (3.27) 

at  5  =  0,  (3.28) 

M  M{s)  as  r  oo,  (3.29) 

where  M{s)  is  the  far-held  prohle  of  M  associated  with  p{s).  Only  one  more  difficulty 
separates  us  from  a  solution  to  this  problem. 

Consider  for  a  moment  the  implications  of  the  material  conservation  of  a*  (3.12)  and 
the  dehnition  of  the  material  conservation  operator  (3.8).  Notice  that  the  future  distri¬ 
bution  of  a*  depends  on  the  radial  advection  by  gravity-inertia  waves.  Our  goal  has  been 
to  neglect  the  transient  adjustment,  but  the  advection  of  a*  by  the  transient  oscillations 
prohibits  this  from  occurring.  This  problem  can  be  circumvented  by  transforming  to  a 
new  coordinate  system  that  makes  u  implicit  in  the  coordinate  transform  and  reduces  the 
material  conservation  to  a  local  conservation. 

3.4»^  Gradient  Balance  Condition  with  :  ential  Radius  Coordinate 

In  order  to  develop  an  alternate  radial  coordinate,  we  must  identify  some  materially 
conserved  property  of  the  how  that  will  allow  us  to  eliminate  u  from  (3.8).  Notice  that 
the  tangential  momentum  equation  (3.4)  may  be  rewritten  in  absolute  angular  momentum 
form 


dM 

ds 


=  Tt 


dM 
'  ds 


=  M 


+  =0.  (3.26) 


^(r»  +  l/r2)=0. 


(3.30) 
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Using  this  conservation  principle,  Schubert  and  Hack  (1983)  defined  the  potential  radius 
coordinate  R  which  is  related  to  physical  radius  by 


=  rz; -h  (3.31) 

R  represents  the  radius  that  a  parcel  in  physical  space  must  be  moved,  while  conserv¬ 
ing  angular  momentum,  for  its  tangential  wind  component  v  to  vanish.  Mathematical 
consistency  requires  that  the  R  coordinate  be  accompanied  by  a  corresponding  S  and  T 
coordinate  in  order  to  distinguish  between  partial  derivatives  at  fixed  r  {d/ds  and  d/dt) 
from  those  at  fixed  R  {d/dS  and  d/dT).  Derivatives  in  either  space  are  related  by 

Applying  (3.32)  to  (3.8)  reveals  that 


-  A 

Dt  ~  dT' 


(3.33) 


Advection  by  u  has  been  eliminated  from  the  calculation  by  making  it  implicit  in  the 
transformation.  The  material  consei'vation  in  physical  space  has  been  converted  into  a 
local  conservation  in  transformed  space  allowing  us  to  rewrite  (3.12)  as 


where  a*  in  the  transformed  coordinate  may  be  written  in  Jacobian  form 


(3.34) 


_  f_  _  dp  ^  d{^r‘^,s)  di^r^,p)  _  di^r^,p) 

C«  "  d{^li^)ds~  di^R:^S)di^r^,s)  ~  di^R‘^,Sy 

With  (3.32)  and  (3.35),  the  invertibility  principle  (3.26)  transforms  to 


(3.35) 


dikr^p) 

di^R^S) 


=  0, 


fv*  = 


dM 

Jr' 


(3.36) 


(3.37) 
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dM 

dS 


(3.38) 


where  M  =  M  +  and  v*  =  {R/r)v  are  the  trainsformed  streamfunction  and  tangential 
velocity,  respectively.  Combining  (3.36)-(3.38),  a  single  expression  for  the  invertibility 
principle  in  terms  of  M  is  obtained 


[/2  _  ^3  A  f  A  (e-S/c.2^) 

[■'  dR  [R^dRjl  ds  \  ds 


f  d'^M  \ 
\dRdS 


) 

)  pcp  \ 


2 

fRdRJ 


=  0. 


(3.39) 


The  upper  and  lower  boundary  conditions  remain  the  same  except  for  being  expressed  in 
transformed  coordinates 


^  =  Tt  at  S  =  St,  (3.40) 

/o  2  dM\  1  /dMV  o  ^  ^  /o 

where  St  =  Cpln(6^T/0o)-  Along  the  outer  boundary  R  =  Rmaxt  a  boundary  condi¬ 
tion  (Dirichlet^)  that  simulates  an  infinite  domain  is  adopted.  If  Rmax  ^  where 
A'  is  associated  with  the  primary  mode  of  the  forcing,  then  a*  and  M,  along  the 
outer  boundary  are  assumed  independent  of  R.  With  this  assumption,  the  invertibil¬ 
ity  principle  at  =  Rmax  is  solved  as  a  one-dimensional  problem  in  5.  Specifically, 
(^""{Rmax^S)  ==  a[Rmax^S)  =  —d'pjdS  is  integrated  to  obtain  p  and  (3.38)  is  integrated 
along  Rmax  to  obtain  Af . 

No  physical  boundary  is  actually  present  at  iZ  =  0;  consequently,  instead  of  defining 
a  boundary  condition  here,  we  will  find  a  limiting  form  of  the  invertibility  principle.  We 
begin  with  the  average  vorticity  within  a  horizontal  circle  of  radius  r 


^Given  a  function  u(a:,y,  z),  the  solution  of  the  Laplacian,  V^u  =  0,  over  a  region  R  with  boundary 
condition  u(x,y,  2)  =  /(a:,y,  2)  on  the  boundary  C  is  called  a  Dirichlet  problem.  If  the  boundary  curve 
C  and  the  boundary  value  function  /  are  reasonably  well  behaved,  then  a  unique  solution  exists  to  the 
problem. 
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2  /■'■ ,  /  w  ^  2u  ^ 

—  (rdr=f+  —  =  f  + 
Jo  r 


2  dM 


imR- 

If  C  approaches  some  finite  limit  Co  as  0,  then  v  =  0{r).  This  being  the  case, 
according  to  the  definition  of  i2,  r  and  R  must  tend  toward  zero  together  (i.e.,  r  =  0{R) 
and  R  =  OCr)).  The  limit  of  (3.42)  as  r  — >•  0  is 


^  2dM\ 


(3.43) 


Recognizing  from  (3.37)  that  lim  dM/dR  =  0,  we  can  apply  L’Hopital’s  Rule^  to 

ft — ^0 

(3.43)  giving 


C  =  /  +  y^  at  R  =  0.  (3.44) 

Substituting  this  expression  into  a  =  iC/f)<^*  along  with  the  definition  of  tr,  the  ideal 
gas  law,  and  the  hydrostatic  relation  (3.38),  we  obtain  the  following  invertibility  relation 
along  R  =  0 


\  2  d^M\  ,  c/c  ^ 


=  0  at  R  =  0. 


(3.45) 


The  upper  boundary  condition  associated  with  this  invertibility  principle  remains  un¬ 
changed 


dM 


=  Tx  at  i?  =  0,  5  =  Sx'y 


(3.46) 


however,  the  lower  boundary  condition  becomes 


M  —  c. 


dM 


=  0  at  iZ  =  0,  5  =  0. 


(3.47) 


(3.45)”(3.47)  represents  the  invertibility  principle  at  i?  =  0. 


\imf{x)/g{x)  results  in  the  indeterminate  form  0/0  or  oo/oo,  then  the  \imf{x)/g{x)  = 
lim/^(a;)/5f'(x)  provided  that  the  latter  limit  exists  (or  is  infinite). 
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3.5  Analyzing  the  Balanced  Fields 


3.5.1  Energetics 


When  developing  a  numerical  model,  it  is  always  prudent  to  analyze  the  energy  of 
the  system  to  ensure  that  it  is  consistent  with  your  formulation  and  the  phenomena  being 
modeled.  The  flux  form  of  the  KE  equation  for  (3.4)-(3.9)  is 


d{aK)  d{ruaK)  dM  ^ 

- 1 - L  fjril -  =  Q 

dt  rdr  dr 


(3.48) 


where  K  =  +  is  the  KE.  Applying  the  definitions  M  =  CpT-^gz  and  a  =  —dp/ds^ 

and  remembering  that  the  flow  is  adiabatic,  it  can  be  shown  that 


diaK)  ,  d[ruaiK  ^  (I))]  ,  d  f  ^dp\  ^ 

-IT* - +  = 


(3.49) 


where  uj  =  Dp/Dt  is  the  vertical  p-velocity.  KE  is  generated  through  thermodynamic 
expansion  and  dissipated  through  work  done  along  the  boundaries  by  pressure.  Similarly, 
the  potential  energy  (PE)  is  given  by 


d{aH)  d{ruaH) 

^  --  aao;  =:  0. 

ot  ror 


(3.50) 


where  H  =  CpT  is  the  enthalpy.  Notice  that  aau  represents  the  conversion  between  KE 
and  PE.  Summing  (3.49)  and  (3.50)  yields  the  total  energy  equation 


d[a{K^H)]  ^\rua{K  -f  M)]  ^ 

dt  *  rdr  d& 


Vs  “  “• 


(3.51) 


To  find  the  total  energy  of  the  system,  we  integrate  over  the  domain  noting  that  at  the 
upper  boundary  dp/dt  =  0  and  at  the  lower  boundary  0  =  0 


:{IC  +  n)  =  {ruaiK  +  M)}^^ 


^ — f'max  ’ 


(3.52) 


where  (}  =  ds, 


/C  =  /  /  aKrdrds 


(3.53) 
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is  the  total  KE  and 


'Ji  =  J J  aHrdrds^ 


(3.54) 


is  the  total  PE.  Notice  that  the  total  energy  is  not  conserved  by  the  system.  During 
the  transient  adjustment,  gravity-inertia  waves  transport  KE  and  PE  across  the  outer 
boundary  r  =  Tmax  and  out  of  the  system.  Using  the  definition  of  a  and  (3.9),  H  may  be 
rewritten 


1  +  kJJ  \pqJ 

The  minimum  H  obtainable  by  adiabatically  redistributing  mass  is 


(3.55) 


PoJ 


(3.56) 


where  p  =  p{s)rdr  is  the  averages  pressure  on  an  isentropic  surface.  The 

difference  T-i  —  H.  is  the  amount  of  PE  that  may  be  converted  into  KE  (i.e.,  the  APE) 


APE 


_  ^OPO  f  f 

~ 1+kJJ 


1+''  f  p\ 


)/ 


i^/^^rdrds. 


(3.57) 


Notice  that  the  APE  is  also  dependent  on  the  outer  boundary.  The  APE,  as  defined  by 
Lorenz  (1955),  is  defined  for  a  finite  horizontally  boundless  domain  (i.e.,  an  atmosphere 
bounding  a  spherical  surface).  As  Vmax  is  increased,  the  APE  will  decrease.  A  similar 
problem  was  encountered  by  Shutts  and  Gray  (1994).  They  found  that  the  outer  bound¬ 
ary  must  be  at  least  a  Rossby  deformation  radius  away  from  the  perturbation.  Model 
sensitivity  studies  show  that  an  outer  boundary  distance  of  2000  km  is  sufficient  for  the 
simulations  to  be  considered  in  the  next  chapter. 

3.5.2  Mass  Removed 

The  intensity  of  the  resulting  balanced  circulations  will  obviously  be  a  function  of  the 
amount  of  mass  that  is  redistributed.  We  could  compare  the  balanced  flow  to  a;  however, 


the  physical  meaning  of  a  is  difficult  to  grasp.  Instead,  we  choose  to  compare  the  balanced 
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flow  to  the  fractional  amount  of  mass  removed  from  a  layer  in  relation  to  the  total  mass 
in  the  layer.  The  amount  of  mass  removed  is  given  by 


Jo  Jo 


(3.58) 


fr  fST  {  2'Kg  ^G^^dr^ds  >  0, 

0  a'  <  0. 

Likewise,  the  total  amount  of  mass  in  each  layer  from  which  mass  was  removed  is  given 

by 


m 


-fl 

Jo  Jo 


ST  (  2'Kg  ^ar^dr^ds  cr'  >  0, 


(3.59) 

^0  JO  [q  ^/  <  0. 

The  ratio  ji  =  m! /rfi  indicates  the  amount  of  mctss  within  a  given  radius  r  used  to  force 

the  balanced  circulation. 


3.6  Applying  the  Appropriate  Numerics 


Using  a  method  similar  to  that  applied  by  Schubert  and  Alworth  (1987),  the  invert- 
ibility  principle  in  the  interior  (3.39)-(3.41)  and  along  i?  =  0  (3.45)-(3.47)  is  discretized 
using  second-order  centered  finite  differences  on  a  uniform  single  64  x  128  grid.  Vertical 
differences  are  applied  at  the  midpoints  of  the  grid  with  ghost  points  outside  the  top  and 
bottom  boundaries.  A  simple  Gauss-Seidel  point  relaxation  scheme  is  used  to  obtain  the 
balanced  solution  except  at  the  outer  boundary.  Since  JA  is  assumed  independent  of  R 
along  the  outer  boundary,  the  invertibility  principle  is  solved  as  a  semi-geostrophic  one¬ 
dimensional  boundary  value  problem.  In  some  instances,  when  the  vertical  or  horizontal 
coupling  is  strong,  a  line  relaxation  method  is  used.  Also,  multigrid  processing,  as  de- 
scribed  by  Pulton  (1989),  is  used  to  accelerate  the  simulations.  The  nonlinearity  of  these 
calculations  is  clearly  evident  in  the  partial  differential  equation  (3.39)  and  lower  bound¬ 
ary  condition  (3.47).  Additional  nonlinearity  is  present  in  (3.39)  due  to  the  dependence 
of  p  on  JA  through  the  hydrostatic  relation.  Despite  the  nonlinearities,  a  unique  solution 
is  attainable. 

Ignoring  the  dependence  of  p  on  A4,  the  discretized  partial  differential  equation  (3.39) 
can  be  expressed  as  a  quadratic.  Schubert  and  Alworth  (1987)  showed  that  if  the  fluid  is 
limited  to  those  cases  where  it  is  statically  and  inertially  stable  (i.e.  cr  >  0  and  (a  >  0) 
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then  a  single  solution  to  the  quadratic  is  obtained.  To  improve  on  this  result  and  improve 
the  model  stability,  the  dependence  o^  p  on  M  is  accounted  for  by  solving  a  Newtonian 
iteration  on  A4.  When  the  mass  perturbation  is  relatively  strong,  the  nonlinearities  may 
cause  the  solution  to  become  unstable.  Under  this  circumstance,  a  continuation  method 
is  employed. 

Using  a  weaker  solution  as  the  initial  approximation  to  the  perturbed  field,  the  relax¬ 
ation  is  continued.  As  many  as  ten  continuations  may  be  used  to  obtain  a  stable  solution. 
Once  a  solution  has  been  obtained,  the  resulting  fields  are  transformed  from  (ii,  S)  to 
(r,  s)  using  cubic  Hermite  splines  with  monotonicity  preserved. 


Chapter  4 


RESULTS 


4.1  Introduction 

In  Chapter  3,  a  nonlinear,  axisymmetric,  gradient  balanced  model  was  developed  to 
compute  the  adjusted  state  of  an  initially  stationary,  continuously  stratified  atmosphere 
following  an  instantaneous  vertical  redistribution  of  mass.  This  chapter  will  present  and 
discuss  results  from  numerous  simulations  conducted  with  the  model.  First,  a  control 
simulation  is  developed  and  compared  to  the  observational  data  presented  in  Chapter  2 
(Section  4.2).  A  series  of  sensitivity  tests  are  then  performed  to  determine  the  effects  of 
variable  forcing,  static  stability,  latitude,  and  cloud-top  cooling  on  vortex  structure  and 
strength  (Section  4.3). 

4.2  Developing  the  Control  Simulation 

In  this  study,  the  control  simulation  serves  a  dual  purpose.  First,  as  with  any  exper¬ 
imental  control,  it  is  used  as  a  base  for  sensitivity  analysis;  however,  for  our  purposes,  it 
will  also  be  used  as  a  model  comparison  to  observational  data.  By  designing  the  control  in 
this  way,  the  physics  of  the  model  can  be  validated  against  observational  evidence  prior  to 
conducting  the  sensitivity  studies.  With  this  in  mind,  the  control  simulation  is  developed 
in  the  following  sections. 

4 ^2 A  Initial  Conditions 

The  basic  state,  initial  conditions  for  the  control  simulation,  as  defined  by  the  tem¬ 
perature  profile  T{s)  (Fig.  3.3a),  are  obtained  from  the  observational  data  in  Chapter  2 
and  are  listed  in  Table  4.1a.  Most  of  the  initial  conditions  are  drawn  from  the  JAN 


46 


sounding  (Fig.  2.7).  Although  this  sounding  is  not  complete,  it  does  show  the  height  of 
the  tropopause  approximately  350  km  from  LCH  and  according  to  BN,  is  representative 
of  other  soundings  far  removed  from  the  convection.  To  simplify  the  analytic  form  of  T{s) 
in  relation  to  the  JAN  sounding,  the  stratospheric  lapse  rate  is  assumed  to  be  isothermal 
throughout  its  depth  while  the  tropospheric  lapse  rate  linearly  increases  (in  isentropic  co¬ 
ordinates)  to  the  surface.  A  surface  temperature  is  chosen  which  reproduces  the  average 
lapse  rate  of  the  LCH  and  SIL  soundings.  Plotted  in  Fig.  4.3  as  the  dashed  curve  is  the 
resulting  T{s)  profile.  This  profile  is  converted  to  pressure  and  then  differentiated  with 
respect  to  entropy  to  obtain  the  initial  a{s)  profile  shown  in  Fig.  4.1a. 


Parameter  Name 

Symbol 

Value 

latitude 

<!> 

30°N 

rotation  rate 

n 

7.292x10-55-1 

top  entropy 

St 

1.0  (100  mb) 

top  temperature 

Tt 

-60C 

tropopause  entropy 

^TR 

0.37  (200  mb) 

tropopause  temperature 

Ttr 

-60C 

tropopause  depth 

^str 

1C 

surface  temperature 

To 

26C 

a 


Parameter  Name 

Symbol 

Value 

perturbation  radius 

r2 

300  km 

smoothing  radius 

n 

240  km 

perturbation  amplitude 

a 

3.6  kgKm-3 

entrainment  amplitude 

e 

0 

perturbation  height 

54 

0.64  (150  mb) 

parameters  from 

53 

0.49  (175  mb) 

Fig. 

52 

0.37  (200  mb) 

5l 

0.09  (700  mb) 

b 

Table  4.1:  Control  settings  of  the  (a)  basic  state  and  (b)  perturbation  parameters.  The 
values  of  entropy  s  =  Cp  ln(0/0o),  have  been  nondimensionalized  by  st  =  Cp  ln(0r/^o)-  To 
further  clarify  the  meaning  of  these  values,  equivalent  pressure  level  have  been  included. 


Similarly,  the  perturbation  initial  conditions,  as  defined  by  the  F{t)  and  G{s)  profiles 
illustrated  in  Fig.  3.4,  are  obtained  from  the  observational  data  in  Chapter  2  and  are 
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listed  in  Table  4.1b.  The  perturbation  radius  r2  wa3  chosen  by  estimating  the  average 
radius  of  the  cloud  shields  observed  by  FM  and  BN.  Both  were  found  to  be  consistent 
with  CLMT’s  composite  values  of  Rc.  Obviously  the  mass  forcing  within  the  outflow  is 
not  constant  across  r2  but  gradually  diminishes  near  the  boundary.  To  account  for  this, 
the  mass  perturbation  is  radially  smoothed  to  zero  over  the  outer  60  km  of  r2.  Since  the 
model  is  impulsively  rather  than  gradually  forced,  the  magnitude  of  the  mass  forcing  a 
is  restricted  to  the  amount  of  mass  available  in  the  basic  state  at  each  level  a(s)  which 
ranges  from  0.112  kgKrn^^  at  the  top  to  3.774  kgKm“^  at  the  surface.  In  the  vertical, 
the  depth  of  the  mass  sink  si  was  chosen  by  considering  both  the  planetary  boundary 
layer  depth  at  LCH  (Fig.  2.7)  and  the  depth  of  surface  convergence  analyzed  by  CLMT 
(Fig.  2.5a).  Likewise,  the  depth  and  vertical  placement  of  the  mass  injection  were  chosen 
by  considering  the  structure  of  the  outflow  at  LCH.  Notice  that  55  is  absent  from  the 
table.  This  parameter  is  only  used  when  cloud-top  cooling  is  included  (i.e.,  e  7^  0).  The 
resulting  vertical  profile  of  G{s)  is  shown  in  Fig.  4.1b.  Both  a{s)  and  i7'(r,  s)  =  aF{r)G{s) 
are  combined  to  produce  cr(r,  5),  the  convectively  perturbed  mass  distribution,  which  for 
an  atmosphere  at  rest  is  simply  (/"“(r,  5)  (Fig.  4.4). 

Notice  that  the  thermodynamic  structure  of  the  troposphere  is  markedly  altered  by 
the  impulsive  redistribution  of  mass.  As  a  clue  to  this  change,  Fig.  4.4  shows  that  the 
center  of  the  a*  anomaly  appears  to  be  below  the  tropopause  (or  the  ba.se  of  the  mass 
injection)  in  pressure  coordinates  but  not  in  isentropic  coordinates.  As  the  atmospheric 
center  of  mass  is  shifted  upwa  rd,  a  hydrostatic  pressure  increase  develops  through  the 
column  extending  from  the  top  of  the  perturbation  to  the  surface  (Fig.  4.2).  This  pressure 
anomaly  not  only  increases  the  tropopause  pressure,  causing  the  injection  to  appear  deeper 
in  isobaric  coordinates,  but  also  adiabatically  warms  the  column,  producing  a  10-15° C 
temperature  change  between  200  and  300  mb  (Fig.  4.3). 

With  the  initial  a*  field  defined  from  a,  the  invertibility  relation  is  solved  and  the 
balanced  state  of  the  control  simulation  obtained  (Figs.  4.5-4.10).  To  explain  the  structure 
of  the  balanced  state,  in  relation  to  the  initial  structure,  we  must  understand,  at  least 
quaditatively,  the  transient  adjustment  that  has  occurred. 
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Figure  4.1:  Control  simulation  initial  (a)  a(5)  and  (b)  G{s)  profiles. 


4.2,2  Balanced  State 

Horizontal  force  imbalances,  produced  by  the  mass  perturbation,  generate  gravity- 
inertia  waves  which  transport  mass  radially.  As  a  result,  cr*  is  materially  advected  outward, 
increasing  the  radial  extent  of  the  anomaly  (compare  Fig.  4.4  with  Fig.  4.5).  Recognizing 
that  r  =  R  for  an  atmosphere  initially  at  rest,  we  can  determine  the  amount  of  displace¬ 
ment  from  the  final  distribution  of  R  (Fig.  4.10).  Compare  these  displacements  with  those 
obtained  by  Shutts  et  al.  (1988)  (Fig.  2.9c).  The  nonlinear  advection  by  gravity-inertia 
waves,  implicit  in  the  coordinate  transformation,  is  clearly  evident.  If  the  displacements 
are  large  relative  to  A',  the  flow  will  adjust  toward  gradient  balance  producing  the  anti- 
cyclonic  anomaly  shown  in  Fig.  4.6a. 

The  outward  transport  of  mass  also  causes  the  initial  hydrostatic  pressure  anomaly 
to  decrease,  adiabatically  cooling  the  troposphere.  However,  it  will  not  cool  back  to  its 
initial  temperature  (Fig.  4.6b)  due  to  the  small  amount  of  mass  that  remains  following 
the  adjustment.  Less  clear  in  Figs.  4.5-4. 6  are  the  vertical  displeicements  (see  Fig.  2.9c) 
produced  during  the  adjustment. 
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Figure  4.4:  Control  simulation  initial  a*  field  in  (a)  isobaric  and  (b)  isentropic  coordinates 
(70  kgKm“^  contour  interval).  Notice  the  change  in  the  perturbation  depth  between  (a) 
and  (b). 


r  (km)  r  (km) 

a  b 

Figure  4.5:  Same  as  Fig.  4.4  except  that  the  fields  are  for  the  adjusted  state.  It  appears 
that  the  perturbation  forces  the  tropopause  downward  along  its  base  and  produces  an 
additional  inversion  along  its  top.  Notice  that  the  sloping  “eye”  wall  structure  on  the 
oncer  edge  of  the  surface  a"'  anomaly  is  similar  to  that  obtained  by  Shutts  et  al.  (1988). 


Figure  4.6:  Control  simulation  adjusted  state  including  the  (a)  tangential  velocity  v 
(1  ms”^  contour  interval)  and  potential  temperature  6  fields  (10  K  contour  interval), 
and  the  (b)  basic  state  (dashed)  and  adjusted  (solid)  temperature  profiles. 
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Figure  4.7:  Control  simulation  (a)  buoyancy  frequency  N  (0.005  contour  interval) 
following  the  adjustment  including  the  (b)  inner  (solid)  and  outer  (dashed)  boundary 
profiles.  Notice  the  sharp  peaks  in  N  above  and  below  the  perturbation. 


Figure  4.8:  Control  simulation  isentropic  relative  vorticity  following  the  adjustment 
(5.0  X  10“^  contour  interval). 
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Figure  4.9:  Control  simulation  aspect  ratio  \' /H'  (contour  interval  is  50). 


r  (km) 

Figure  4.10:  Control  simulation  potential  radius  R  following  the  adjustment  (100  km 
contour  interval).  Contours  represent  surfaces  of  constant  angular  momentum.  Notice 
that  the  upper  displacements  are  much  greater  than  those  at  the  surface. 


54 


Unbalanced  vertical  pressure  gradients  produced  by  the  mass  perturbation  generate 
vertically-propagating  gravity  waves.  It  would  appear  from  the  governing  equations  (3.4)- 
(3.8)  that  vertical  advection  is  not  present  in  the  model;  however,  recognizing  that 


we  see  that  vertical  advection  is  implicit  in  the  material  derivative.  Thus,  the  adiabatic 
ascent  and  descent  of  vertically  propagating  gravity  waves  produces  a  warm-cold  couplet 
(Fig.  4.6b).  Amazingly,  the  vertically  and  horizontally  propagating  gravity  waves  work  in 
harmony  to  produce  a  temperature  and  velocity  distribution  that  maintains  thermal  wind 
balance. 

Horizontal  temperature  gradients,  essential  to  the  maintenance  of  thermal  wind  bal¬ 
ance,  result  from  the  bowing  of  isentropic  surfaces.  Apparent  in  the  isentropes  of  Fig.  4.6a 
is  a  lens  structure  similar  to  that  obtained  by  Gill  (1981)  and  Shutts  et  al.  (1988)  (Figs.  2.8 
and  2.9b).  In  terms  of  energy,  the  lens  represent  the  fraction  of  APEq  remaining  in  the 
system  following  the  adjustment  or  in  other  words,  the  balanced  APE.  Notice  that  the 
changes  in  the  temperature  distribution  also  alter  the  static  stability  (Fig.  4.7). 

Increased  vertical  temperature  gradients  along  the  horizontal  boundaries  of  the 
anomaly  result  in  increased  static  stability  while  the  decreased  temperature  gradient 
within  the  interior  of  the  anomaly  results  in  decreased  static  stability.  Laxge  static  stability 
gradients  may  cause  gravity  wave  trapping  or  ducting.  Although  the  transient  adjustment 
is  not  solved  explicitly  in  this  model,  one  might  speculate  that  as  the  peaks  in  N  de¬ 
velop,  the  propagation  of  gravity-inertia  waves  (and  the  adjustment)  will  be  increasingly 
restricted  to  the  horizontal.  A  comparison  of  Figs.  4.6a  and  4.7a  shows  that  the  adjust¬ 
ment  spreads  to  greater  depths  as  the  vertical  gradient  of  N  weakens.  Restricting  the 
adjustment  depth  causes  intense  shear  zones  to  develop  above  and  below  the  anomaly. 

Strong  shear  is  also  produced  in  the  horizontal  as  shown  in  the  distribution  of  isen¬ 
tropic  relative  vorticity  C  (Fig.  4.8).  A  strong  gradient  of  cyclonic  vorticity  is  found  along 
the  outer  edge  of  the  anomaly.  In  contrast,  the  interior  is  dominated  by  nearly  constant 
anticyclonic  vorticity.  This  structure  agrees  well  with  both  Gill’s  (1981)  and  Shutts  et 
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al.’s  (1988)  results  (Figs.  2.8  and  2.9).  Along  the  outer  edge  of  the  perturbation,  the  flow 
is  becoming  more  inertially  stable  while  the  interior  is  becoming  less  inertially  stable.  Ra¬ 
dial  displacements  are  inhibited  by  the  inertial  stability  such  that  the  largest  displacements 
or  strongest  adjustment  is  experienced  by  parcels  within  the  interior  of  the  perturbation 
(Fig.  4.10). 

Recall  that  Gill  (1981)  and  Shutts  et  al.’s  (1988)  linear  formulations  showed  that  the 
aspect  ratio  of  the  lens  is  given  by  A/.ff  =  N/f.  In  regions  were  N  varies  significantly 
from  that  of  the  basic  state  or  C  becomes  large  in  comparison  to  /,  the  aspect  ratio  is  more 
accurately  represented  by  the  rotational  form  A'/ =  N(r,  s) / ^ {(  pjg^  4  9 

shows  that  most  of  the  domain,  except  for  the  interior  of  the  anomaly,  is  dominated  hy  X/H 
which  for  the  stratosphere  is  approximately  200.  Decreased  inertial  stability  within  the 
interior  of  the  perturbation  causes  a  four  fold  increase  in  the  aspect  ratio  {X'/H^  «  800). 
Thus,  for  a  given  value  of  A'  increases  within  the  interior  resulting  in  a  shift  of  the 
adjustment  from  the  velocity  to  the  mass  field  (Schubert  et  ah,  1980). 

At  the  surface,  gravity-inertia  waves  transport  meiss  inward  to  fill  the  void  left  by 
the  mass  sink;  however,  Fig.  4.10  indicates  that  lateral  displacements  are  minimal  (in 
comparison  to  the  upper  perturbation).  As  pointed  out  by  Shutts  et  al.  (1988),  radial 
inflow  is  limited  by  the  conservation  of  angular  momentum.  If  the  flow  was  allowed  to 
converge  onto  the  axis,  conservation  of  angular  momentum  would  produce  an  infinite 
tangential  velocity;  consequently,  to  fill  the  sink,  mass  must  also  descend  (Fig.  2.9c). 
Adiabatically  descending  gravity  waves  produce  the  warm  core  low  (Fig.  4.6a  and  b)  and 
sloping  “eye  wall”  frontal  structure  (Fig.  4.5)  observed  by  Shutts  et  al.  (1988). 

4  >2.3  Comparison  to  Observations 

At  this  point,  we  can  compare  the  results  of  the  model  to  the  observations  presented 
in  Chapter  2  and  discuss  the  validity  of  the  model  physics.  First,  consider  the  similarities 
between  the  model-derived  wind  field  (Fig.  4.6a)  and  the  perturbation  wind  fields  obtained 
by  FM  (Fig.  2.1-Fig.  2.4).  The  simulated  anticyclone  heis  a  minimum  velocity  of -13.5  ms“^ 
at  500  km  which  is  approximately  a  third  less  than  the  analyzed  20  ms"^.  Two  factors 
contribute  to  this  difference.  First,  the  latitudes  of  the  systems  studied  by  FM  were  5^-10® 
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farther  north  than  the  30°  used  in  the  simulation.  As  will  be  shown  in  a  later  section, 
increasing  latitude  increases  the  strength  of  the  balanced  circulation.  Second,  the  amount 
of  forcing  may  be  stronger  in  the  observational  studies.  Recall  that  the  model’s  forcing  is 
limited  by  the  amount  of  mass  initially  available  in  the  basic  state.  If  the  model  were  time 
dependent,  mass  could  be  supplied  continuously  resulting  in  greater  forcing  and  stronger 
balanced  circulations.  At  the  surface,  a  warm  core  cyclone  develops  with  a  maximum 
velocity  of  11.5  ms“^  at  156  km.  CLMT’s  analysis  does  show  the  development  of  cyclonic 
vorticity  at  the  surface  prior  to  and  during  the  systems  formation;  however,  the  mesoscale 
cyclone  which  develops  in  the  model  has  more  in  common  with  oceanic  systems  such 
as  polar  lows  (Rasmussen,  1985;  Shapiro  et  al.,  1987)  and  tropical  cyclones  (Shea  and 
Gray,  1973)  than  it  does  with  those  systems  observed  over  land.  More  detailed  simulation 
including  the  effects  of  downdrafts  are  needed  to  fully  understand  the  mid-tropospheric 
and  surface  features  analyzed  by  CLMT. 

Interesting  similarities  between  the  model  results  and  observations  can  also  be  seen 
in  the  adjusted  temperature  profile  (Fig.  4.6b),  and  the  soundings  at  LCH  and  JAN 
(Fig.  2.7).  The  warm-cold  couplet  produced  by  the  lens  has  a  vertical  structure  that  is 
nearly  identical  to  the  temperature  distributions  in  the  soundings.  Model  warm  and  cold 
anomalies  are  5.3°C  and  4.2°C,  respectively,  as  compared  to  the  2-3°C  and  13°C  of  the 
soundings.  Recall  that  FM  observed  a  1“2°C  warm  anomaly  and  a  6-8°C  cold  anomaly. 
Evidently,  the  model  warming  is  too  strong  while  the  model  cooling  is  much  too  weak. 
The  discrepancies  can  be  partially  accounted  for  as  follows. 

Recall  from  Fig.  4.3  that  the  nonhydrostatic  pressure  increase,  caused  by  the  initial 
mass  perturbation,  generates  a  deep  layer  of  warm  air,  A  remnant  of  this  effect  is  present 
in  the  warm  anomaly.  Since  the  actual  transient  vertical  motions  which  generate  the  mass 
perturbation  are  nonhydrostatic,  it  is  believed  that  this  hydrostatic  effect  is  an  artifact  of 
the  model  and  not  an  actual  physical  process;  consequently,  the  magnitude  of  the  warm 
pool  is  overestimated  by  approximately  1-2°C. 

Discrepancies  in  the  cold  anomaly,  on  the  other  hand,  are  caused  by  the  lack  of  addi¬ 
tional  cooling  atop  the  perturbation  produced  by  entrainment  and  radiative  processes  not 
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included  in  the  model.  The  impact  of  these  processes  will  be  considered  in  later  sensitivity 
studies.  Entrainment  and  turbulent  mixing  also  impact  the  temperature  gradient  within 
the  injection. 

Notice  that  the  observed  ;  liabatic  gradient  within  the  anvil  at  LCH  is  not  repro¬ 
duced.  Since  the  model  flow  is  aydrostatic,  turbulent  motions  are  not  generated;  thus,  the 
observed  homogeneous  layer  is  not  reproduced.  The  model  does,  however,  form  the  sharp 
inversions  that  are  observed  above  and  below  the  injection. 

Interestingly,  the  double  inversion  structure  is  dynamically  linked  to  the  a*  field. 
Inversions  or  gradients  in  static  stability  are  associated  with  gradients  in  a*.  Fig.  4.5 
shows  that  the  inversion  splits  at  the  perturbation.  One  branch  represents  the  depressed 
tropopause  while  the  other  is  formed  by  the  cold  pool  atop  the  lens.  Several  sensitiv¬ 
ity  studies  will  now  be  examined  to  determine  the  impact  of  varying  the  perturbation 
parameters  in  Table  4.1.  The  forcing  parameters  a  and  T2  will  be  considered  first. 

4.3  Comparing  the  Control  to  Sensitivity  Simulations 

^,3,1  Forcing  Changes 

Fig.  4.11  shows  the  relationship  between  a  and  the  percentage  of  mass  removed  and 
injected  within  two  perturbation  radii  r2.  In  the  forcing  diagrams  that  follow,  a  will  be 
shown  on  the  bottom  of  the  figures  and  the  percentage  of  mass  removed  along  the  top. 
In  addition,  the  control  simulation  is  identified  by  the  blackened  data  points.  Fig.  4.12 
illustrates  the  effect  of  varying  a  on  the  velocity  extremes  and  the  parcel  displacements 
necessary  to  produce  those  extremes. 

As  one  might  expect,  increases  in  forcing  produce  stronger  vortices;  however,  not 
as  obvious  are  the  nonlinear  growth  rates.  When  a  =  0.5,  the  anticyclone  is  four  times 
stronger  than  the  cyclone;  whereas,  when  a  =  3.6,  the  anticyclone  strength  has  decreased 
to  1.2  that  of  the  cyclone.  It  appears  that  the  anticyclone  growth  approaches  a  linear  limit 
while  the  cyclone  growth  continues  to  accelerate.  This  is  again  consistent  with  angular 
momentum  conservation.  When  d  «  4,  the  magnitudes  of  the  vortices  are  identical; 
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Figure  4.11:  Relationship  between  the  perturbation  amplitude  a  and  the  percentage  of 
mass  removed  from  a  cylinder  of  radius  2r2,  bounded  by  the  surface  and  5i.  The  blackened 
data  points  indicate  the  values  used  in  the  control  simulation. 

however,  the  parcel  displacements  are  not.  The  displacement  of  the  velocity  maodmum  is 
only  about  half  that  of  the  minimum. 

Likewise,  forcing  increases  produce  more  energetic  vortices  (Fig.  4.13a).  Notice  that 
the  APEo  is  an  order  of  magnitude  greater  than  the  balanced  APE  and  KE,  indicating 
significant  gravity-inertia  wave  generation.  When  the  forcing  is  weak,  about  5.8%  of 
the  APEo  is  retained  and  another  9%  is  converted  into  balanced  KE  (Fig.  4.13b).  The 
remaining  85.2%  radiates  beyond  the  domain  as  gravity-inertia  waves  during  the  transient 
adjustment.  Apparently,  nonlinearities  also  affect  the  partitioning  of  the  energetics.  As  the 
forcing  is  increased,  the  fraction  of  the  APE  remaining  decreases  to  4,2%  while  the  fraction 
converted  to  balanced  KE  decreases  to  8%.  Forcing  increases  decrease  the  inertial  stability 
of  the  system  and  reduce  the  efficiency  of  the  adjustment  such  that  less  disturbance  energy 
is  partitioned  to  the  balanced  state  of  the  system. 

Increases  in  r2  also  impact  the  amount  of  mass  forcing.  Variations  in  r2  are  equivcilent 
to  changes  in  the  size  of  the  convective  system.  Below  about  200  km  radius,  the  cyclone  is 
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Figure  4.12:  Variation  of  the  (a)  velocity  extremes  and  the  (b)  radial  displacements  needed 
to  produce  the  velocity  extremes  in  relation  to  the  forcing.  The  blackened  data  points  are 
the  same  as  in  Fig.  4.11. 
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Figure  4.13:  Relationship  between  the  energetics  and  the  forcing  including  the  (a)  inte¬ 
grated  APEo,  APE,  and  KE;  and  the  (b)  percentage  of  APEq  partitioned  to  the  balanced 
APE  and  KE.  The  blackened  data  points  are  the  same  as  in  Fig.  4.11. 


stronger  than  the  anticyclone  (Fig.  4.14a);  whereas,  larger  scale  systems  produce  anticy¬ 
clones  that  are  stronger.  At  a  radius  of  600  km,  the  magnitude  of  the  anticyclone  is  nearly 
a  third  greater  than  that  of  the  cyclone.  As  the  systems  become  smaller  than  100  km,  the 
magnitude  of  the  vortices  asymptotically  decreases  to  zero. 


0  200  400  600  800  1000  0  200  400  600  800  1000 
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a  b 

Figure  4.14:  Variation  of  the  (a)  velocity  extremes  and  the  (b)  percentage  of  APEq  par¬ 
titioned  to  the  balanced  APE  and  KE  as  a  function  of  the  perturbation  radius  r2.  The 
blackened  data  points  are  the  same  as  in  Fig.  4.11. 

As  systems  grow,  more  of  the  APEq  is  converted  into  balanced  energy  (Fig.  4.14b). 
A  doubling  of  the  radius  from  200  to  400  km,  increases  the  fractional  APE  from  2%  to  6% 
and  the  fractional  KE  from  4%  to  7%.  According  to  Schubert  et  al.  (1980),  as  increases 
with  respect  to  A,  less  of  the  APEq  is  partitioned  to  gravity-inertia  waves  and  more  to 
the  balanced  flow.  In  other  words,  the  adjustment  continually  shifts  from  the  mass  field 
to  the  velocity  field.  The  next  section  considers  the  impact  of  changing  static  stability  on 
the  structure  of  the  adjusted  state. 

4^3.2  Static  Stability  Changes 

If  the  updraft  fails  to  penetrate  into  the  stratosphere,  the  outflow  will  come  to  rest 
below  the  tropopause.  Fig.  4.15  shows  the  adjusted  state  of  a  simulation  that  is  identical 
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to  the  control  except  that  the  50  mb  perturbation  is  now  located  immediately  beneath 
the  tropopause. 

Overall,  the  balanced  state  of  the  system  weakens  relative  to  the  control.  The  anti¬ 
cyclone  magnitude  has  decreased  by  27%  to  -9.87  ms“^  with  an  80  km  decrease  in  radial 
extent  to  420  km.  Also,  there  is  a  3°C  reduction  in  the  cold  anomaly  and  a  2.5^C  reduction 
in  the  warm  anomaly.  In  contrast  to  the  control,  the  tropopause  is  forced  upward  with  no 
secondary  temperature  minimum.  By  maintaining  a  50  mb  perturbation  depth,  the  mass 
forcing  is  identical  to  the  control  simulation;  therefore,  the  variations  can  only  be  due  to 
the  change  in  static  stability  incurred  by  moving  the  injection  from  the  stratosphere  to 
the  troposphere. 

Variations  in  static  stability  are  measured  in  terms  of  the  Brunt  Vaisala  or  buoyancy 
frequency  N  =  yj pg^/acp.  Fig.  4.16b  indicates  that  the  basic  state  N  is  nearly  constant 
in  both  the  troposphere  and  stratosphere  with  a  large  discontinuity  at  the  tropopause. 
The  phase  speed  of  an  internal  gravity  wave  within  a  continuously  stratified  fluid  is 
proportional  to  AT.  Thus,  as  the  static  stability  decreases,  gravity-inertia  waves  propagate 
more  slowly  and  the  adjustment  occurs  over  shorter  horizontal  scales  (i.e.,  A  =  Cn/ f 
decreases;  compare  Figs.  4.10  and  4.19).  In  addition,  according  to  angular  momentum 
conservation,  the  balanced  flow  will  be  weaker  since  the  horizontal  displacements  are 
smaller.  Recall  that  the  aspect  ratio  is  also  proportional  to  N , 

In  this  case,  A'/ff'  is  approximately  200  which  is  about  double  that  of  the  troposphere 
(A/jEf  «  100).  Thus,  for  a  given  H\  the  decreased  N  within  the  troposphere  results  in  a 
A'  which  is  smaller  than  that  of  the  control  simulation.  Large  values  of  X^/H^  atop  the 
anomaly  show  the  influence  of  increased  static  stability  above  the  tropopause.  Under¬ 
standing  the  impact  of  N  on  the  adjustment  is  further  clarified  by  considering  the  changes 
in  the  balanced  state  that  result  from  gradually  moving  the  perturbation  through  the 
tropopause. 

As  the  perturbation  comes  into  contact  with  the  tropopause  and  the  increased  static 
stability  of  the  stratosphere  above,  the  anticyclone  magnitude  rapidly  increases  from  9  to 
15  ms^^  (Fig.  4.20a).  For  vertical  separations  greater  than  if,  between  the  perturbation 
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Figure  4.16:  Same  as  Fig.  4.7  except  that  the  50  mb  perturbation  is  located  immediately 
beneath  the  tropopause. 


Figure  4.17:  Same  as  Fig.  4.8  except  that  the  50  mb  perturbation  is  located  immediately 
beneath  the  tropopause. 


Figure  4.18:  Same  as  Fig.  4.9  except  that  the  50  mb  perturbation  is  located  immediately 
beneath  the  tropopause. 


Figure  4.19:  Same  as  Fig.  4.10  except  that  the  50  mb  perturbation  is  located  immediately 
beneath  the  tropopause. 


Figure  4.21:  Variation  of  the  (a)  minimum  velocity  and  the  (b)  percentage  of  APEq 
partitioned  to  the  balanced  APE  and  KE  as  a  function  of  the  average  stratospheric  static 
stability  iV.  The  blackened  data  points  are  the  same  as  in  Fig.  4.11. 
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and  the  tropopause,  the  magnitude  of  the  anticyclone  as  a  function  of  altitude  becomes 
a  constant.  An  alternate  method  to  test  the  influence  of  N  on  the  adjustment  involves 
maintaining  the  vertical  position  of  the  perturbation  while  varying  the  basic  state  N.  For 
a  doubling  of  iV,  the  magnitude  of  the  anticyclone  increases  by  about  4  ms“^  (Fig.  4.21a). 

Notice  that,  as  the  anticyclone  intensifies,  the  partitioning  of  APEq  to  the  balanced 
APE  and  KE  decreases  (Figs.  4.20b  and  4.21b).  This  again  shows  that  the  partitioning 
of  the  energy  is  related  to  the  inertial  stability  of  the  system.  Recall  that  A  is  not  only 
dependent  on  N  but  on  /  as  well.  In  the  next  section,  the  influence  of  changing  latitude 
on  the  adju: :  ment  is  examined. 

^.5.5  Latitude  Changes 

Fig.  4.22  shows  the  adjusted  state  of  a  simulation  identical  to  the  control  except  that 
the  perturbation  has  been  displaced  south  to  10°N.  Again,  as  in  the  previous  simulation, 
the  adjusted  state  is  weaker  than  the  control.  The  cold  and  warm  temperature  anomalies 
have  decreased  by  5.5  and  3°C,  respectively,  but  with  the  same  overall  vertical  structure. 
In  addition,  the  magnitude  of  the  anticyclone  has  diminished  by  over  a  half  to  -6.63  ms”^; 
however,  in  contrast  to  the  simulation  testing  variability  due  to  changing  iV,  the  radius  of 
the  minimum  velocity  has  increased  to  615  km  (Fig.  4.26).  Figs.  4.22-4.26  clearly  show 
that  the  horizontal  scale  A  of  the  adjusted  state  has  increased  while  the  vertical  scale  H 
has  decreased. 

In  this  case,  the  aspect  ratio  of  the  basic  state  stratosphere  has  doubled  in  comparison 
to  the  control  {\/H  ^  400).  More  striking,  however,  is  the  tenfold  increase  within  the 
perturbation  {X^/H^  ^  4000).  At  low  latitudes,  where  /  is  small.  A'  is  extremely  sensitive 
to  changes  in  the  inertial  stability  of  the  vortex  (Fig.  4.24).  Thus,  as  the  inertial  stability 
decreases  within  the  interior  of  the  anomaly,  A^  increases  dramatically  causing  the  adjusted 
flow  to  weaken  (Fig.  4.22).  In  addition,  the  decreased  inertial  stability  allows  for  greater 
horizontal  displacements  of  adjusting  parcels  (Fig.  4.26). 

Several  sensitivity  simulations  have  been  conducted  to  study  the  impact  of  changing  / 
on  the  kinematics  and  energetics  of  the  system.  Fig.  4.27a  shows  that  the  influence  of  /  is 
greatest  where  f3  =  df  /dcj)  is  large.  A  20°  change  in  latitude  within  the  subtropics  results 
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Figure  4.27:  Variation  of  the  (a)  velocity  extremes  and  the  (b)  percentage  of  APEq  par¬ 
titioned  to  the  balanced  APE  and  KE  as  a  function  of  the  latitude.  The  curves  do  not 
extend  to  the  equator  due  to  the  limitations  of  the  /-plane  approximation.  Note  that 
APEq  is  not  a  function  of  latitude,  but  has  a  constant  value  1.23xlO^^J.  The  blackened 
data  points  are  the  same  as  in  Fig.  4.11. 

in  a  7  ms“^  change  in  the  anticyclone  and  a  3  ms*”^  change  in  the  cyclone.  Interestingly, 
33  the  equator  is  approached,  the  cyclone  becomes  stronger  than  the  anticyclone.  As 
the  dependence  of  A'  on  (  increases  at  low  latitudes,  the  inertially  stable  cyclone  will 
receive  more  balanced  energy  than  the  less  inertially  stable  anticyclone.  Schubert  and 
Alworth  (1987)  showed  that  cyclone  strength  is  further  enhanced  if  the  mass  source  occurs 
over  a  deeper  layer. 

For  convective  systems  near  the  equator,  where  /  approaches  zero,  values  of  A  become 
extremely  large;  thus,  the  energy  partitioning  shifts  almost  entirely  from  the  balanced  flow 
to  gravity  wave  generation  (Fig.  4.27b).  A  change  in  latitude  from  30°N  to  lO^'N  results 
in  an  amazing  69%  and  80%  decrease  in  the  KE  and  APE,  respectively.  Thus,  for  tropical 
convective  systems  to  develop  intense  mesoscale  anticyclones  aloft,  they  must  be  much 
larger  than  their  mid-latitude  counterparts.  Earlier  in  the  chapter,  it  was  pointed  out 
that  the  cold  anomaly  is  much  too  weak.  We  return  to  this  subject  and  consider  the 
moderating  effects  of  entrainment  and  radiation. 
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4.4  Considering  the  Effects  of  Entrainment  and  Radiation 


Figure  4.28:  Profile  of  G(s)  including  an  additional  mass  sink  along  the  top  of  the  outflow. 

Probably  the  most  significant  discrepancy  between  the  model  results  and  the  obser¬ 
vational  data  is  the  large  difference  in  the  magnitude  of  the  cold  anomaly.  It  appears 
that  one  or  more  additional  diabatic  processes  are  affecting  the  cloud-top  temperature. 
The  most  likely  candidate  would  be  longwave  radiative  cooling;  however,  as  pointed  out 
in  Chapter  2,  the  downward  flux  of  heat  from  the  stratosphere  by  entrainment  is  also  im¬ 
portant.  For  example,  Shutts  and  Gray’s  (1994)  plume  model  produces  a  small  region  of 
cloud-top  cooling  in  the  absence  of  radiative  cooling.  They  suggest  that  the  “existence  of 
such  a  cool  layer  is  consistent  with  a  small  net  mass  sink  created  by  the  entrainment  of  air 
into  the  convective  plume  during  its  overshoot  phase.”  Recall  from  the  model  formulation 
that  provisions  were  made  in  the  vertical  perturbation  profile  G{s)  to  add  an  additional 
mass  sink  at  the  top  of  the  injection. 

Such  a  sink  will,  in  a  manner  similar  to  the  method  used  by  Fritsch  and  Brown  (1982), 
reproduce  the  diabatic  cooling  that  occurs  at  cloud  top.  By  choosing  e  ^  0,a,  small  amount 
of  mass  is  removed  from  the  layer  above  and  drawn  down  into  the  injection.  Choosing 
e  =  0.1  produces  the  maximum  possible  perturbation.  A  perturbation  top  of  S5  =  0.64 
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(125  mb)  yields  an  entrainment  depth  of  1.4  km  which  is  consistent  with  the  penetration 
depth  of  an  entraining,  overshooting  cloud-top.  The  modified  perturbation  profile  is  shown 
in  Fig.  4.28. 

Keep  in  mind  that  the  perturbation  depth  is  decreased  in  comparison  to  the  control 
profile,  causing  the  mass  forcing  to  increase  from  67%  to  85%.  Thus,  cloud-top  diabatic 
effects  not  only  affect  the  cold  anomaly  but  the  gradient  adjustment  process  as  well.  The 
adjusted  state  of  this  simulation  is  shown  in  Fig.  4.29. 


e 

Figure  4.30:  Variation  of  the  cold  pool  temperature  as  the  magnitude  of  the  cloud-top 
mass  sink  e  is  increased.  The  blackened  data  point  are  the  same  as  in  Fig.  4.11. 

The  decreased  depth  of  the  perturbation  along  with  the  increased  mass  forcing  has 
strengthened  the  anticyclone,  in  comparison  to  the  control,  by  0.8  ms“^.  Note  also  that 
a  shallow  cyclonic  lens,  with  a  magnitude  less  than  2  ms~^,  forms  above  the  anticyclone. 
Shutts  and  Gray’s  (1994)  plume  simulations  produce  a  similar  vortex.  Cyclonic  motion, 
associated  with  the  vertical  advection  of  cyclonic  momentum,  is  observed  near  the  center 
of  hurricane  outflow;  however,  similar  features  have  not  been  observed  in  the  outflow 
of  MCSs.  This  simulation  suggests  that  cyclonic  vorticity  may  be  generated  within  the 
outflow  through  diabatic  cooling  and  gradient  adjustment.  Keep  in  mind  that  the  forcing 
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used  to  simulate  the  cooling  is  extremely  simplified.  In  reality,  the  radius  of  the  entraining 
updraft  is  much  less  than  r2  while  the  radiative  cooling  occurs  over  a  depth  much  shallower 
than  55  5/^  • 

Fig.  4.30  shows  the  magnitude  of  the  cold  pool  as  a  function  of  e.  The  cold  pool 
temperature  does  decrease  by  an  additional  2.5°C;  however,  the  resulting  temperature 
anomaly  remains  much  less  than  that  observed  by  BN,  suggesting  that  even  stronger 
forcing  is  present.  No  change  is  observed  in  the  magnitude  of  the  warm  anomaly. 


Chapter  5 


SUMMARY  AND  CONCLUSIONS 


5.1  Summary 

After  examining  the  observational  data  of  mesoscale  convection  presented  in  Chap¬ 
ter  2,  the  assertion  was  made  that  the  anticyclonic  outflow  within  the  upper-troposphere 
and  lower-stratosphere  can  be  partially  explained  by  the  gradient  adjustment  of  the  dis¬ 
turbed  mass  field.  To  verify  this  hypothesis,  a  nonlinear  balanced  model  was  developed 
in  Chapter  3  from  the  axisymmetric,  quasi-static  primitive  equations  to  examine  the  gra¬ 
dient  balanced  state  of  a  hypothetical  convective  system.  Specifically,  we  idealized  the 
convection  as  an  impulsive  transfer  of  mass  from  the  lower-troposphere  to  a  layer  near 
the  tropopause  within  a  continuously  stratified  environment.  Using  isentropic  and  trans¬ 
formed  radius  (Schubert  and  Hack,  1983)  coordinates,  we  obtained  an  analytic  solution 
similar  to  the  classic  Rossby  invertibility  problem.  Numerical  simulations  were  conducted 
in  Chapter  4  to  verify  the  model’s  physics  and  to  test  the  model’s  sensitivity  to  the  mass 
forcing,  static  stability,  uititude,  and  cloud-top  diabatic  forcing.  Several  conclusions  can 
be  drawn  from  the  results  of  the  numerical  simulations. 

5.2  Conclusions 

1)  In  general,  the  results  of  the  control  simulation  suggest  that  the  formation  of  the 
observed  anticyclone  can  be  explained  by  the  gradient  adjustment  of  the  convectively 
disturbed  mass  field.  The  model  produces  an  anticyclone  with  accompanying  thermal 
structure  that  is  similar  in  magnitude  and  form  to  that  observed  by  FM  and  BN. 

2)  Other  diabatic  effects  alter  the  adjustment.  When  cloud-top  cooling  (i.e.,  en¬ 
trainment,  radiation,  etc.)  is  included,  the  magnitudes  of  both  the  anticyclone  and  cold 
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temperature  anomaly  are  increased.  In  addition,  a  weak  cyclone  forms  in  the  region  of 
cooling.  Although  this  feature  hcts  not  yet  been  observed  it  does  appear  in  more  complex 
numerical  simulations  (Shutts  and  Gray,  1994). 

3)  Lower-stratospheric  injections  produce  a  kinematic  and  thermodynamic  structure 
that  agrees  with  theory  (Danielsen,  1982)  and  is  similar  to  that  observed  (FM;  Danielsen, 
1982;  and  BN).  In  contrast,  when  the  injection  level  is  within  the  upper-troposphere,  the 
magnitudes  of  the  anticyclone  and  temperature  anomalies  are  much  weaker.  In  addition, 
no  secondary  temperature  minimum  is  produced  at  the  base  of  the  outflow.  Previous 
models  have  neglected  the  impact  of  the  stratosphere  resulting  in  the  formation  of  anti¬ 
cyclones  that  are  weaker  and  have  less  horizontal  extent  (Shutts  et  al.  1988;  Shutts  and 
Gray,  1994).  The  simulations  presented  here  clearly  show  the  importance  of  including  a 
stratosphere  in  any  study  of  deep,  penetrative  convection. 

4)  Large  increases  of  static  stability  along  the  boundaries  of  the  outflow  may  result  in 
gravity  wave  ducting.  Numerical  simulations  of  a  two-dimensional  squall  line  conducted 
by  Tripoli  and  Cotton  (1989a,  b)  show  that  gravity  wave  trapping  along  the  base  of  the 
outflow  was  enhanced  when  surface  long-wave  heating  was  included.  It’s  possible  that  the 
heating  enhances  the  inversion  at  the  base  of  the  outflow  and  increases  the  static  stability. 

5)  Most  of  the  initial  available  potential  energy  is  partitioned  to  gravity-inertia  oscil¬ 
lations.  Only  a  fraction  is  partitioned  to  the  balanced  KE  of  the  adjusted  state.  Schubert 
et  al.  (1980)  demonstrated  that  the  energy  efficiency  is  dependent  on  A,  showing  that 
“for  small-scale  initial  disturbances  in  the  pressure  field  most  of  the  initial  energy  ends 
up  in  gravity-inertia  waves.”  The  amount  of  partitioning  is  still  a  matter  of  contention. 
Schubert  et  al.’s  (1980)  results  show  that  it  is  less  than  10%  while  Shutts  and  Gray  (1994) 
have  obtained  values  greater  than  30%.  The  results  of  this  study  show  that  it  is  close  to 
8%.  Future  research  is  needed  in  this  area  to  determine  what  individual  factors  influence 
the  partitioning. 

6)  The  efficiency  of  balanced  energy  generation  decreases  as  the  anticyclone  inten¬ 
sifies.  Using  an  Eliassen  balanced  vortex  model  to  study  tropical  cyclones,  Hack  and 
Schubert  (1986)  found  that  the  generation  of  balanced  KE  is  a  function  of  the  inertial 
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stability  of  the  vortex.  Intensification  of  the  anticyclone  leads  to  reduced  inertial  stability 
and  decreased  balanced  energy  generation.  Increased  forcing  would  eventually  cause  the 
cyclone  magnitude  to  surpass  that  of  the  anticyclone,  resulting  in  the  increased  inertial 
stability  of  the  system  and  greater  efiiciency  of  balanced  energy  generation. 

5.3  Future  Research 

It  is  important  to  point  out  that  this  model  represents  a  first  approximation  to  the 
dynamic  processes  occurring  within  mesoscale  convective  outfiows.  Many  important  phys¬ 
ical  processes  have  been  neglected  in  this  formulation,  not  because  they  are  unimportant, 
but  because  they  limit  our  ability  to  find  analytic  solutions  to  the  adjustment.  Advance¬ 
ments  in  our  understanding  of  the  adjustment  problem  will  come  from  increasingly  more 
complex  models  which  include  the  effects  neglected  here. 

For  instance,  a  clear  limitation  of  this  model  is  its  lack  of  moisture.  Moist  processes 
must  be  included  in  future  models  to  account  for  changes  in  the  static  stability  (Lalas, 
D.  P.,  and  F.  Einaudi,  1974;  Durran  and  Klemp,  1982)  and  the  temperature  distribution 
within  the  outflow.  The  temperature  profile  is  also  affected  by  cloud-top  entrainment  and 
turbulent  mixing  within  the  layer.  Furthermore,  turbulence  results  in  eddy  diffusion  which 
affects  the  shape  of  the  perturbation  lens  (Gill  et  al.’s,  1979;  Hedstrom  and  Armi,  1988). 
To  simulate  the  turbulence,  future  models  must  be  three-dimensional  and  nonhydrostatic. 
Differential  radiative  heating  across  the  cloud  layer  further  enhances  turbulent  mixing, 
requiring  a  radiative  transfer  model  with  accurate  representations  of  cloud  liquid  and  ice 
water  contents  (Starr  and  Cox,  1985a  and  b;  Ackerman  et  al.  1988).  Increased  under¬ 
standing  of  the  energy  partitioning  dictates  that  future  models  account  for  the  downscale 
cascade  of  energy  and  the  transition  to  quasi-two  dimensional  turbulence  (Lilly,  1988). 
To  understand  the  instabilities  present  in  these  flows,  asymmetries  (Hedstrom  and  Armi, 
1988)  and  vertical  shear  must  be  added  (Shutts,  1987,  1994).  Not  only  can  an  increased 
understanding  of  the  adjustment  process  be  gained  from  improved  numerical  models  but 
from  more  detailed  observational  studies  as  well. 

Doppler  radar,  wind  profilers,  and  improved  satellite  imagery  lend  themselves  to 
modern  studies  of  the  anticyclonic  anomaly.  This  should  include  analysis  similar  to  FM’s 
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to  see  what  impact  mesoscale  convection  has  on  current  large-scaJe  NWP  models.  Recall 
that  in  the  introduction  we  asked  the  question:  What  errors  do  mesoscale  convective 
outflows  in  the  upper-troposphere  and  lower-stratosphere  introduce  into  our  large-scale 
NWP  model  forecasts? 

Certainly,  as  this  study  shows,  increasing  the  magnitude  of  the  system  will  increase 
the  error.  However,  for  systems  of  similar  size,  the  error  introduced  by  the  anticyclonic 
anomaly  will  be  greatest  at  higher  latitudes  for  those  systems  strong  enough  to  inject  mass 
into  the  stable  layers  of  the  lower  stratosphere.  Clearly  the  results  of  this  and  past  studies 
indicate  the  need  to  increase  our  understanding  of  this  process  in  order  to  improve  future 
convective  parameterization  schemes. 
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